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Abstract: We present a complete study of AS = 2 and AB = 2 processes in the left- 
right model (LRM) based on the weak gauge group SU(2)l x SU(2)r x U(1)b-l- This 
includes £k, AMk, AM q , ^4g L , AT q with q = d,s and the mixing induced CP asymmetries 
Si[)K s an d S^ff,. Compared to the Standard Model (SM) these observables are affected 
by tree level contributions from heavy neutral Higgs particles (H°) as well as new box 
diagrams with Wr gauge boson and charged Higgs (H^) exchanges. We also analyse the 
B — > A S) d7 decays that receive important new contributions from the Wl — Wr mixing 
and exchanges. Compared to the existing literature the novel feature of our analysis 
is the search for correlations between various observables that could help us to distinguish 
this model from other extensions of the SM and to obtain an insight into the structure of 
the mixing matrix V R that governs right-handed currents. Moreover, we perform the full 
phenomenology including both gauge boson and Higgs boson contributions. We find that 
even for M H o ~ M H ± ~ 0(20) TeV, the tree level H° contributions to AF = 2 observables 
are by far dominant and the contributions to B — >• X q j can be very important, even 
dominant for certain parameters of the model. While in a large fraction of the parameter 
space this model has to struggle with the experimental constraint from ek, we demonstrate 
that there exist regions in parameter space which satisfy all existing AF = 2, B — > Xs^J, 
tree level decays and electroweak precision constraints for scales Myy R ~ 2 — 3 TeV in the 
reach of the LHC. We also show that the S^k s - tension present in the SM can be 
removed in the LRM. Simultaneously Br(i? —> X s ^f) can be brought closer to the data. 
However, we point out that with the increased lower bound on M\\r R , the LRM cannot 
help in explaining the difference between the inclusive and exclusive determinations of 
\V U b\, when all constraints are taken into account, unless allowing for large fine-tuning. 
Finally we present a rather complete list of Feynman rules involving quarks, gauge bosons 
and Higgs particles. 
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1 Introduction 

A very important property of the Standard Model (SM) regarding flavour violating pro- 
cesses is the left-handed structure of the charged current interactions reflecting the maximal 
violation of parity observed in low energy processes. Left-handed charged currents encode 
at the level of the Lagrangian the full information about flavour mixing and CP violation 
represented compactly by the CKM matrix. Due to the GIM [1] mechanism this structure 
has automatically profound implications for the pattern of FCNC processes that seems to 
be in remarkable accordance with the present data within theoretical and experimental un- 
certainties [2, 3], bearing in mind certain anomalies [2, 4-12], in particular in CP- violating 
observables which are discussed below. 

As the SM is expected to be only the low-energy limit of a more fundamental theory it is 
conceivable that at very short distance scales parity could be a good symmetry implying the 
existence of right-handed (RH) charged currents. Prominent examples of such fundamental 
theories are left-right (LR) symmetric models on which a rich literature exists. Indeed, 
LR symmetric models were born 37 years ago [13-17]. Early papers mainly cover the 
examinations of two special cases, known as "manifest" scenario [17] and "pseudo-manifest" 
scenario [18-21], which are characterised by no spontaneous and fully spontaneous CP 
violation, respectively. The right-handed counterpart of the CKM matrix appears then in 
a special form being either identical to or the complex-conjugate of the CKM matrix up 
to certain phases as e.g. summarised in [22]. The phenomenology of both scenarios has 
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widely been studied in the literature [23-26]. By now in both scenarios strong constraints 
on the heavy charged gauge boson mass Mw R ^ 4TeV have been obtained [23, 27, 28] 
from the constraints on the Kl — K$ mass difference, CP violation in kaon decays and the 
neutron electric dipole moment, making these scenarios difficult to access in direct searches 
at the LHC. In addition the "pseudo-manifest" scenario has been ruled out by both the 
appearance of light Higgs triplets [29] and the correlation of Ek and sin(2/3) [30]. This 
means that the right-handed mixing matrix must be different from the CKM matrix in 
order to reach agreement with experiment. Motivated by this fact more general studies on 
CP violation have been performed in [22, 31, 32]. More recent extensive analyses of many 
observables in the LR symmetric framework can be found in e.g. [28, 33-36]. 

Theoretical interest in models with an underlying SU(2)l x SU(2)r global symmetry 
has also been motivated by Higgsless models [37-40]. Moreover, the recent phenomeno- 
logical interest in having another look at the right-handed currents in general originated 
from tensions between inclusive and exclusive determinations of the elements of the CKM 
matrix \V u b\ and |V^,|. As pointed out and analysed recently in particular in [41-43], 
the presence of right-handed currents could either remove or significantly weaken some of 
these tensions, especially in the case of \V u b\. The implications of these findings for many 
observables within an effective theory approach have been studied in [44]. 

Yet an effective theory approach, as interesting as it may be, involves a number of 
unknown couplings that limit the predictive power of the theory and in particular does 
not allow to correlate low-energy high precision observables to high energy processes be- 
ing already explored in a new domain of energy at the LHC. In this context we refer to 
[45, 46] where extensive analyses of LR symmetric models have been performed for early 
LHC data. Therefore it is of interest as a preparation for new discoveries both through 
high energy processes and high precision experiments in this decade to perform a detailed 
phenomenological analysis in a concrete class of models with right-handed currents, in 
particular models with LR symmetry. As manifest and pseudo-manifest LR models have 
already been ruled out we term the model with extended gauge group but without exact 
P or C symmetry Left-Right Model (LRM). 

The goal of the present paper is to analyse the well measured FCNC observables related 
to the particle-antiparticle mixings K — K° and B® s — B® g in this NP scenario. For a RH 
scale in the reach of the LHC the off-diagonal mixing amplitudes M[ 2 (i = K, d, s) receive 
important and often dangerous tree level contributions from neutral heavy scalar particles 
(H°) present in this model. Additionally box diagrams with exchanges of a heavy W'^ and 
charged Higgs (-ff ) exchanges and right-handed couplings of the light are present. 
Similarly the Wl — Wr mixing can have a significant impact on the B — > decay that 
often puts a severe constraint on extensions of the SM. Also heavy charged Higgs (H^) 
exchanges can contribute and in fact these contributions cannot be neglected as often done 
in the literature. 

We would like to know whether this class of very interesting models can be made 
consistent with all existing data for RH scales as low as M\y R ~ (2 — 3) TeV, which is still 
consistent with direct collider searches, while solving various anomalies observed in the 
quark sector. 
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As there have been other analyses of particle-antiparticle mixing [22, 23, 27, 28, 31, 
33, 47-53] and B — > X s ^y [54-62] in LR models, it is mandatory for us to state what is new 
in our paper: 

• First of all we perform a simultaneous analysis of the most interesting AF = 2 
observables in the K and Bd s meson systems in conjunction with the decays B — > X s ^f 
and B — > X^- The analysis is performed in a general framework of the LRM and 
hence without making from the beginning particular assumptions on the structure of 
the right-handed mixing matrix or equality of gauge couplings. That is, in contrast 
to previous analyses we do not assume a certain specific form for by restricting its 
parameters, but search for its structure by using the bounds from tree-level decays, 
AF = 2 observables, B — > X q ^ decays and imposing constraints from electroweak 
precision tests. This strategy differs from the existing literature which dominantly 
considered bounds on the parameters of LR models, in particular on the masses of 
W r± and H°. 

• In this manner we are led to new structures of the V R matrix that are still rather 
simple and allow to monitor transparently which anomalies observed in the quark 
sector can be solved in these models. In this context we search for correlations 
between various observables. 

• We perform the full phenomenology including both gauge boson and Higgs boson 
contributions finding that even for M H o ~ M H ± ~ O(20) TeV, the tree level H° 
contributions to AF = 2 observables are by far dominant and the contributions 
to B — > Xgj can be very important and even dominant for certain parameters of the 
model. In this context we include the known QCD corrections. 

• We analyse the issue of the element \V u b\ in this specific model with right-handed 
currents, pointing out that with the increased lower bound on Mw R , these models 
cannot help in explaining the difference between the inclusive and exclusive determi- 
nations of \ V u b\, when all constraints are taken into account, unless allowing for large 
fine-tuning of parameters. 

• We investigate a soft lower bound on the heavy Higgs mass. 

• We present a collection of Feynman rules necessary for the analysis of all flavour 
violating processes, in particular AF = 2 transitions considered in the present paper. 
These rules could also be useful for collider physics. 

Our paper is organised as follows. In section 2 we summarise briefly the main in- 
gredients of the LRM. In section 3 we present the effective Hamiltonians for K° — K°, 
B^ — B^ and B® — B® mixings and we calculate the most interesting observables such as 
the CP-violating parameter £k, the mass differences AMk and AM q , the CP-asymmetries 
Ag L (q = d, s), S^k s an d and the width differences Ar g . In section 4 we present the 
analysis of the B — > X s ^ decays including CP-violating asymmetries. In section 5 we face 
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tree level decays including B + — > t + u t . In section 6 we summarise the existing constraints 
on the LRM from the electroweak precision tests. In section 7 we outline our strategy for 
the numerical analysis. In this context we review the existing anomalies in the flavour data 
and present a simple analytical expression for the right-handed matrix V R that allows to 
see how these anomalies can be solved in a correlated manner for different values of \V u b\. 
In section 8 we present the bounds on the electroweak sector of the LRM from electroweak 
precision tests. In section 9 a general study of the right-handed mixing matrix is pre- 
sented. Subsequently in section 10 a detailed numerical analysis of particle-antiparticle 
mixing observables and of B — > X S) dfy decays including tree-level constraints is performed. 
In section 11 we derive a soft lower limit on the heavy Higgs masses. In section 12 a 
brief comparison of the LRM with other models is presented. We summarise our results in 
section 13. In the appendices we provide a more detailed description of the LRM and the 
symmetry breaking mechanism, the Higgs sector and numerical insights into the structure 
of LR contributions to AF = 2 operators. Furthermore we provide an extensive list of 
Feynman Rules for the LRM. 

2 Models with Left-Right gauge symmetry 

In this section we give a brief description of the LRM. We restrict our presentation to the 
key properties, which allows us to set our notation. More details on LR symmetric and 
asymmetric models can be found e.g. in [28, 31, 34] and references therein. 

2.1 Gauge group and fermion content 

Among the most popular new physics models are LR extensions 

517(3) C x SU{2) L x SU(2) R x U(1) B -l (2.1) 

of the SM gauge group, as they allow for a restoration of parity symmetry at high energies. 
Note that the spontaneous breaking of parity does not have to be connected to the break- 
down to the SM gauge group and can take place at some much higher scale. Therefore we 
do not restrict ourselves to the study of the parity (or alternatively charge conjugation) 
symmetric case but consider the generic case with independent gauge couplings gh,9R- 

The left-handed fermions are embedded as SU (2) l doublets and SU (2)_r singlets, while 
the right-handed fermions are SU{2)l singlets and SU(2)r doublets: 

Ql 1 : I ~ [■■■ ■ l . i . Qh 1 "," I -(3,1,2,-^ , (2.2) 





3, 

L L =[ ^ I -(1,2,1,-1), L R = j -(1,1,2,-1) , (2.3) 

with the quantum numbers given in brackets corresponding to the gauge group in (2.1). 
We see that instead of the seemingly arbitrary U(1)y charges in the SM, the fermionic U (1) 
charges are now given by their B — L quantum numbers. The resulting electric charges are 

Q = T 3L + T 3R + . (2.4) 
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2.2 Higgs sector and spontaneous symmetry breaking 

The spontaneous symmetry breaking of LR models takes place in two steps. 

Step 1: At a high scale k r ~ O(TeV) SU(2) R x U(1) B -L is broken to the SM 
hypercharge gauge group U(1)y- 

SU(2) R xU(l) B -L^U(l) Y . (2.5) 

The details of this breaking are model-dependent. The two simplest possibilities introduce 
either two scalar doublets or triplets, however also more complicated SU(2)l^ r representa- 
tions are phenomenologically viable. It turns out that quark flavour phenomenology does 
not depend on the particular structure of the Higgs sector (for more details see appendix 
C). Hence from now on we concentrate on the triplet model, which is appealing in the 
neutrino sector as it naturally generates TeV scale Majorana masses for the right-handed 
neutrinos [63-66] through the VEV of A R defined below. Consequently the light neutrino 
masses are suppressed by the TeV scale see-saw mechanism. 

In the triplet model, the symmetry breaking SU(2) R x U(1)b-l U(l)y is achieved 
by a Higgs triplet A R 



which develops a VEV 



<A R > -I 00 .). (2.7, 



The most recent experimental direct Wr searches find roughly M\\r R > 1.5 — 2 TeV and 
this can only be satisfied with k r much larger than the EWSB scale v. This hierarchy 
of scales implies that all LR effects can be expanded in powers of the small dimensionless 
parameter 

e = v/ KR . (2.8) 

Throughout our phenomenological analysis we keep contributions of up to C(e 2 ), which 
constitute the leading corrections relative to the SM result. 

Keeping up the possibility of straightforwardly incorporating the limits of manifest P 
or C symmetry we also introduce an SU (2) l triplet 

^=( 5+L 'f (1,3,1,2), (2.9) 



51 -5+/V2 



(Al) = _ i6 n , (2.10) 



whose VEV 



is constrained to be kl <^ O(eV) in order not to generate large Majorana masses for the 
left-handed neutrinos [63-65]. 
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Step 2: The second step of symmetry breaking is then achieved by the bidoublet 

-(1,2,2,0), (2.11) 




whose vacuum expectation value 




breaks SU(2)l x SU(2)r to its diagonal subgroup SU(2)y: 

SU(2) L x SU(2) R -> SU(2) V . (2.13) 
This step of spontaneous symmetry breaking takes place at the scale 

v = + K /2 _ 174 QeV . (2.14) 

Together with the spontaneous breaking SU(2) R x U(1)b-l ~^ U(l)y in (2.5) at the scale 
kr 3> v the VEV of <fi leads to the standard electroweak symmetry breaking (EWSB) 
SU(2) L x U(1) Y -> £/(1)q so that finally 

Si7(2) L x ST/(2) B x *7(1) B _ L -> U(l) Q . (2.15) 

As we see later in (2.20) in order to obtain a mass splitting between up and down type 
quarks k ^ k 1 is required. More explicitly requiring the hierarchy m& <C mt to be natural 
implies «' <C k. However, we do not assume that k/k = mt/nrib as done in some papers. 
In fact such large values of k/k' are disfavoured by electroweak precision observables, in 
particular by A b FB (see section 8 for details). In our analysis we confine the values to 
1 < k/k' < 10. It should be stressed again that the limit k = k' is not allowed, which can 
be seen explicitly from the divergent behaviour of several observables. 

For our choice of scalar fields the Higgs Lagrangian is then given by 

£ffi g gs = Tr[(^A L )t(^A L )] + Tr[( J D At A /? )t(^ A/? )] + Tr[(D^(D^)] + 1/(0, A L , A R ) , 

(2.16) 

where the covariant derivatives are 

= cV + ig L (^ Llt ■ f)<P - ig R <p(^Rn ■ r) , (2.17) 

D V A (L,R) = ^M A (i,fl) + W(L,R) W (L,R)/i ' T, &(L,R) + W B^( LtR ) . 

and the Higgs potential, as being used in our analysis, is given in appendix A. 
2.3 Gauge sector after electroweak symmetry breaking 

After performing the two steps of symmetry breaking the gauge boson mass matrices can 
be constructed and diagonalised. We summarise our results in appendix B. In the process 
of electroweak symmetry breaking only the gluons and the photon remain massless, while 
the and the Z boson and their heavy counterparts W'^ and Z' acquire masses as 
given in appendix B. The spontaneous symmetry breaking additionally introduces mixing 
between the light and heavy gauge bosons. In the case of the light SM like bosons 
this introduces right-handed couplings at C(e 2 ). The Z boson couplings are also modified 
with respect to the SM but they do not enter the present analysis. 
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2.4 Yukawa interactions and fermion masses 

The most general renormalisable Yukawa coupling of the quark fields with our choice of 
Higgs fields is given by 

£y U k = -yijQu&QRj - yijQu^Qm + n - c - > ( 2 - 18 ) 

where eft = 020* °"2 an d flavour indices i,j = 1, 2, 3. Note that the quantum numbers of the 
other scalar fields in the theory preclude their direct coupling to quarks. 1 The resulting 
fermion mass matrices read 

{M u )ij = v{Y u ) ij , {M d )ij = v{Y d )ij , (2.19) 

where 

(Y u )ij = VijC + Vijse~ ta , {Y d )ij = yijse ia + y^c , (2.20) 

and s = k' /v and c = n/v. These matrices are diagonalised by the bi- unitary transforma- 
tions 

M*°* = U ] L M U U R , (2.21) 
M d diag = D\M d D R , (2.22) 

where Ul r, Dl^r are unitary matrices connecting the flavour and mass eigenstates of 
quarks. 

As we discuss in more details below, the extended Higgs sector with respect to the 
SM leads to flavour changing neutral Higgs couplings already at the tree level. While the 
flavour violating couplings of the light SM-like Higgs are highly suppressed and therefore 
irrelevant for the study of AF = 2 observables and K and B decays, the new heavy Higgses 
lead to dangerously large effects in FCNC observables. In particular the structure of the 
Yukawa coupling of the Higgs bidoublet <p m (2.18) leads to couplings of the down-type 
quarks proportional to the up-type quark mass matrix. Since up and down masses are not 
diagonalised simultaneously, this leads to flavour changing couplings already at the tree 
level, see table 18 for details. 

2.5 Parameter counting 

Having introduced the LR model, let us now count the parameters present in the theory. 
The gauge sector is parametrised by the gauge couplings 

9s, 9L, 9r, 9 , (2.23) 

i. e. one additional parameter relative to the SM. 

The Higgs potential, see appendix A, introduces several new parameters. However we 
see below that for our phenomenological considerations effectively only four parameters 
appear to be relevant. Setting kl = these are 

v = + K /2 j s = k'/v, k r , M b , (2.24) 

1 In the triplet model Al,r couple to the left- and right-handed leptons respectively, generating Majorana 
mass terms M VL R ~ k,l,r- 
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where the first three parametrise the Higgs VEVs 2 and Mh is at leading order the common 
mass of the heavy Higgses H®, H® and (for more details see appendix C). Note that 
that the SM VEV v is given in terms of n and k! . 

In the most general case, the Yukawa couplings yij and y^j are arbitrary complex 
matrices, i. e. contain each 9 real parameters and 9 phases. However, not all of these 
parameters are physical but some can be removed by unitary transformations under the 
flavour symmetry SU(3)q l x SU(3)q r . Finally we are left with the six quark masses and 
two mixing matrices in the LH and RH sectors respectively: 

V L = U{D L , V R = U R D R . (2.25) 

Adopting the standard CKM phase convention, where the 5 relative phases of the quark 
fields are adjusted to remove 5 complex phases from the CKM matrix V L , we have no more 
freedom to remove the 6 complex phases from V R . In the standard CKM basis V R can be 
parametrised as follows [44] 

V R = D v V R tf D , (2.26) 

where V R is a "CKM-like" mixing matrix, containing only three real mixing angles and one 
non-trivial phase. The diagonal matrices Djjd contain the remaining CP-violating phases. 
Choosing the standard parametrisation for V R we have 

( C12C13 S12C13 Sl3e _i ^ 

-S12C23 - Cl2S 2 3Sl3e trt> Cl 2 C 2 3 ~ Sl2S23Sl3e^ S23C13 I , (2.27) 



yR 



\ S12S23 - Cl2C23Sl3e^ -S23C12 ~ Su^Sl^ C23C13 



and 

Djj = diag(l, e^, e ^3 ) , D D = diag(e^, e^, . (2.28) 



3 AF = 2 transitions 
3.1 Preliminaries 

In what follows we use conventions and notation from our papers on various extensions 
of the SM. An easy comparison with the results for AF = 2 observables in the SM, 
the Littlest Higgs model with T-parity (LHT) [67], the Randall-Sundrum scenario with 
custodial protection (RSc) [68] and the SM4 [69] is facilitated in this manner. 

In the LR models the effective Hamiltonian for AF = 2 observables is constructed by 
evaluating three classes of diagrams: 

• The standard box diagrams with quarks and WlWl, WrWr and WlWr exchanges 3 . 
Among the NP contributions involving Wr only the latter matter. 

2 The phase a in the Higgs potential appears in the analytic expressions below. However we eventually 
set it to zero as the factor exp(ia) always multiplies V R and cancels out in all expressions for FCNC 
processes. On a more technical note: the factor exp(ia) always multiplies V R , which is a unitary matrix 
with six phases. Therefore we are always able to absorb a through a redefinition of all phases. Recall that 
a unitary matrix cannot have more than six independent phases. 

3 In what follows in order to make the expressions more transparent it is useful to denote W and W' by 
Wl by Wr, respectively, even if they differ by 0(e 2 ) corrections. 
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• Box diagrams with charged Higgs and gauge boson exchanges. As in LR models 

have masses in the multi TeV range, among NP contributions involving only 
the ones with and Wl matter. 

• Tree level neutral heavy Higgs exchanges. These contributions are problematic unless 
the masses of new neutral Higgs particles are significantly larger than the one of Wr. 
In first approximation the masses of are equal to the masses of the neutral Higgs 
bosons in question. 

While in the SM only one operator contributes to each AF = 2 transition, in the model 
in question there are 8 such operators of dimension 6. Consequently the renormalisation 
group (RG) QCD analysis becomes more involved and due to the LR structure of the new 
operators QCD corrections play a much more important role in new physics contributions 
than in the SM contributions. 

In what follows, after listing all contributing operators we summarise the effective 
Hamiltonian for AF = 2 transitions. To this end we give the formulae for the Wilson 
coefficients at the matching scales and we summarise the RG QCD corrections and the 
results for the hadronic matrix elements. Subsequently we give the final formulae for 
the basic mixing amplitudes in terms of Wilson coefficients at the high scale and the 
effective parameters P{ that encode perturbative and non-perturbative QCD effects. Then 
we discuss in detail the general and special anatomy of LR contributions and compare the 
operator structure to the one found in other NP scenarios. We end this section with listing 
the relevant observables. 

3.2 Local operators 

The contributing operators can be split into 5 separate sectors, according to the chirality of 
the quark fields they contain. For definiteness, we shall first consider operators responsible 
for the K°—K° mixing. The operators belonging to the first three sectors (VLL, LR and 
SLL) read [70] : 

Q\ R (K) = (s a lfi P L d a )(sP^P R dP), 
Q\ K {K) = (s a P L d a )(s?PRdP), 

Qf LL (K) = (s a P L d a )(sPp L d?), 

Ql hh (K) = (s a <y^P L d a )(^a^P L d^) , (3.1) 

where = |[7/n,7i/] and Pl,r = 2(1 =F 75)5), and summation over the colour indices 
a,f3 = 1,2,3 is understood. The operators belonging to the two remaining sectors (VRR 
and SRR) are obtained from Q^ LL and Qf LL by interchanging Pi and Pr. In the SM only 
the operator QY LL (K) is present. The operators relevant for B q (q = d, s) are obtained by 
replacing in (3.1) s by b and d by q. 
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3.3 Effective Hamiltonian 

The effective Hamiltonian for AF = 2 transitions can be written in a general form as 
follows 

r 2 m 2 

,AF=2 U F 1V1 W L 



cff 



An 2 



(3.2) 



where Qi are the operators given in (3.1) and Ci{jj) their Wilson coefficients evaluated 
at a scale fi which we specify below. In what follows we collect the Wilson coefficients 
of these operators separating the contributions from box diagrams with Wl and Wr ex- 
changes from box diagram charged Higgs contributions and tree level neutral Higgs 
H® 2 contributions so that 



Ci = AboxCi + A H ±Ci + A H oCi . 



(3.3) 



These coefficients depend sensitively on the elements of the matrices V L and V R through 
[49] 

\? B (K) = V£*V£ , \f B (B q ) = Vi*V* , (3.4) 
where A, B = L, R, q = d, s and i = u,c, t. 

3.3.1 Wilson coefficients from gauge boson box diagrams 



u, c, t 



W 



± 




W± 




it, c, t 



it, c, t 




W 




u, c, t 



u, c, t 




lu, c, t 




Figure 1. Feynman diagrams for contributing gauge boson box diagrams 

Calculating the diagrams in figure 1, we find the following contributions to the Wil- 
son coefficients relevant for the K° — K° system at the relevant matching scales = 
0{M w ,m t ) and fi R = 0{M Wr ) 



i,j=c,t 

A Box C^ R , K) = >H R (K)Xf L (K)S LR (x u Xj , (3) , 

i,j=u,c,t 

A Box CY RR (v R ,K)= £ Af R (^)Af R (^)5 RR (x,,x,), 

i,j=c,t 



where we introduced the ratios 

2 



X j 



\M Wl 



111; 



Xj 



{m Wr J ' 



( Sw V 

\cwsrJ 



(3.5) 
(3.6) 
(3.7) 

(3.8) 



- 10 - 



Note that at 0(e 2 ) there are no corrections to the Wilson coefficient of the SM Q7 LL 
operator. In the LR symmetric limit gi = gR the factor r reduces to r = 1. 
The loop functions are given as follows 

Shh(xi,Xj) = F(xi,Xj) + F(x u ,x u ) - F(xi,x u ) - F(xj,x u ) , (3.9) 

S LK (xi,Xj,(3) = 2f3r^/XiXj [(4 + x i x j (3)I 1 (x i , xj, (3) - (1 + p)I 2 (x i} Xj, P)] , (3.10) 

SxR(xi,Xj) = f3r 2 S LL (xi,Xj) , (3.11) 

F(xi,Xj) = - [(4: + XiXj)l2(xi,Xj,l) - SxiXjh (xi,Xj,l)\ , (3-12) 



with 

h(xu Xj ,P) ^ -I- (/ -n- '/,, • (3.13) 



Xj ln(xj) ^ln(jg) 

(1 - Xi )(l - XiP)( Xi - Xj) +[Z ^ J) (1 - p)(\ - Xi p)(l - Xj p) 



Tt a\ x 2 ln(xi) ln(/3) 

h(xi, Xj ,p) = 1 ' , + — . 3.14 

(1 - Xi)(l - XiP){Xi - Xj) (1 - P){1 - XiP)(l - XjP) 

The remaining coefficients vanish in the absence of QCD corrections but as we discuss 
below they are generated by QCD effects. In obtaining the results in (3.5) and (3.7) we 
have used the unitarity of the matrices V L and V R or equivalently the GIM mechanism to 
eliminate the A^ L and terms. The GIM mechanism does not apply to the case of LR 
contributions. In the case of B® — B® mixing we just have to replace K by B q . 

The results given above were obtained by calculating all box diagram contributions 
in 't Hooft-Feynman gauge keeping both gauge boson and Goldstone boson contributions. 
We confirm the results in the literature [28, 49, 50, 71, 72]. 

While 5ll and Srr are gauge independent, this is not the case of S'lr as pointed 
out in [71]. As anticipated in that paper and explicitly demonstrated in [50, 72] a gauge 
independent result for Slr is obtained by including vertex and self-energy corrections to 
the tree-level physical Higgs H® and H® exchanges that we discuss subsequently. The 
diagrams relevant for the cancellation of the gauge dependence of 5lr are the ones that 
include the vertices HfG + G ~ but for consistency all vertex and self-energy corrections 
should be included. Detailed analyses in [50, 72] show that the main role of these addi- 
tional contributions is the restoration of the gauge invariance of S'lr without any relevant 
modification of the 't Hooft-Feynman gauge result given above. Consequently we neglect 
these contributions in our analysis. 

3.3.2 Wilson coefficients from charged Higgs box diagrams 

Among various box diagrams with charged Higgs exchanges the one involving and 
W L (see left hand side in figure 2) and the corresponding diagram with Wl replaced by 
the Goldstone bosons are by far the most important ones. The remaining box diagrams 
involving only or and can be safely neglected. A similar comment applies to 
box diagrams with heavy neutral Higgs particles. Calculating the relevant diagram shown 
in figure 2 yields 

A H± C^ H ,K)= J2 ^ R (K)Xf L (K)S? R (x h Xj,p H ). (3.15) 

i,j=u,c,t 
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li, c, t ' F 



i 1 1A, C, t 





Figure 2. Feynman diagrams for Higgs contributions 



with the master function S^i defined as 



S^ R (xi,Xj,(3 H ) = 2u(s)Ph\I Xi{nH)xj{iiH) [xiXjIi(xi,Xj, fi H ) ~ h(xi,Xj, Ph)] • (3.16) 
Here the function u{s) and j3n are given as follows 



u(s) 



1 — 2s 2 



(3.17) 



We have confirmed the result of [49] except that in that paper s <C 1 has been chosen. Our 
choice of the matching scale is explained below. Note that the quark masses in the first 
factor in (3.16) have to be evaluated at /j,h as this factor arises from the Yukawa couplings 
of . The remaining masses have to be evaluated at scales as discussed in the context of 
QCD corrections below. In table 6 we provide values of the masses at different high scales. 

3.3.3 Wilson coefficients from tree level Higgs exchange 

Calculating the tree diagrams with neutral Higgs exchanges (right hand side in figure 2), 
we find that the contributions of the Higgs bosons Hi and fl§ to Cf\ L and Cf^ R cancel 
each other in the limit of M H o = M H o = Mjj. Consequently in the case of non-degenerate 
Higgs masses these contributions are C(e 4 ) and can be neglected. On the other hand in 
the same limit taking into account the overall factor in (3.2) we find 



A H oC 2 LR (^,i^) 



16tH 



V2M 2 H G F 



u{s) ^ R (K)Xf J (K) y /x i (fj. H )x j ( t i H ), (3.18) 



i,j=u,c,t 



with other Wilson coefficients vanishing at = 0{Mh) in the absence of QCD corrections. 
The quark masses have to be evaluated at hh- This result agrees with [28] except that in 
that paper s <C 1 has been chosen. As seen in (3.17) e. g. for s = 0.5 the additional factor 
u(s) provides an enhancement of a factor of 4. In the case of B® — B® mixing K should be 
replaced by B q . A H [)Cj LR is generated by QCD effects as discussed below. 

3.4 QCD corrections and hadronic matrix elements 

The complete analysis of AF = 2 processes requires the inclusion of the QCD renormali- 
sation group evolution. The local operators have to be evolved from their respective high 
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scales down low energy scales at which the hadronic matrix elements are evaluated by 
lattice methods. A complication arises in the model in question as three rather different 
high scales [iw ^ M-R *C hh ar e involved. Before addressing this problem let us recall a 
very efficient method [73, 74] for the inclusion of all these QCD effects in the presence of 
a single high scale which we denote by /Xj n . 

Instead of evaluating the hadronic matrix elements at the low energy scale we can 
choose to evaluate them at the high scale /i; n at which heavy particles are integrated out. 
Thus the amplitude for M — M mixing (M = K,Bd,B s ) is given simply by 

G 2 M 2 

A(M -+ M) = (M\n A J= 2 \M) = £ Cf(fMn) (M\Qi(Wn)\M) . (3.19) 

Here the sum runs over all the operators listed in (3.1). The matrix elements for B^ — Bd 
mixing are for instance given as follows [73, 74] 

(B° d \QU^)\B° d ) = 2 -m 2 Bd F 2 Bd P?(B d ) , (3.20) 

where the coefficients P?(Bd) collect compactly all RG effects from scales below Hi Q as well 
as hadronic matrix elements obtained by lattice methods at low energy scales. Analytic 
formulae for these coefficients are given in [73] while recent applications of this method 
can be found in [44, 75, 76]. As the Wilson coefficients Cj(//j n ) depend directly on the 
loop functions, tree diagram results and fundamental parameters of a given theory, this 
formulation is very transparent and interesting short distance NP effects are not hidden 
by complicated QCD effects. The numerical values for the coefficients Pf(B q ) and P^{K) 
that we require for our analysis are given below. 

The question then arises how to generalise this method to the case at hand which 
involves three rather different high scales. There are three types of contributions for which 
the relevant high energy scales attributed to the coefficients quoted above differ from each 
other: 

• The SM box diagrams involving Wl and the SM quarks. Here the scale is chosen to 
be 0(mt) as in [73]. 

• Tree diagrams mediated by neutral heavy Higgs exchanges. In this case we take 

= 15 TeV as the initial scale for the RG evolution. 

• The only problematic cases are the contributions from Wr and that appear in 
box diagrams together with much lighter Wl and the SM quarks. Here the correct 
procedure would be to first integrate out Wr and . Subsequently one would 
construct an effective field theory not involving them as dynamical degrees of freedom. 
We believe that in view of several unknown parameters in the LR models such a 
complicated analysis would be premature. Therefore we choose [Ir as the matching 
scale for box contributions involving Wr. For diagrams involving Higgs particles we 
set the high scale to be As the dominant effects from the included RG evolution 
stem from scales below fa, this procedure should sufficiently well approximate the 
true result. 
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Now let us turn to the question of scales in the quark masses. In box contributions we 
use rrii(mi) for i = c , b , t and mj(2 GeV) for light quarks. An exception are the masses in 
the overall factor in contribution in (3.16), as discussed previously, which originate in 
the Yukawa couplings of quarks to . Here similarly to the tree level exchange of heavy 
neutral Higgs particles (3.18) quark masses should be evaluated at mj. 

Having the initial conditions for Wilson coefficients at a given high scale /ii n we can 
calculate the relevant M — M amplitude by means of (3.19) provided also the corresponding 
hadronic matrix elements are known at this scale. As seen in (3.20) these matrix elements 
are directly given in terms of the parameters P"(if), P^(Bd) and P^(B S ). Explicit ex- 
pressions for the latter in terms of RG QCD factors and the non-perturbative parameters 
Bf(fiL) are given in equations (7.28)-(7.34) in [73] with \xl denoting the low energy scale 
to be specified below. 

The parameters Bf(fi£,) are subject to considerable uncertainties. They can be ex- 
tracted from the results of [77, 78]. The parameters B^hl) (i = 1, . . . , 5) are given in the 
basis used by Ciuchini et al [79]. Both papers provide the values of the parameters B{ in 
the NDR scheme used in [73]. The conversion to our operator basis (3.1) is given by 

BY LL M = BY RR = B 1 ^ L ), (3.21) 

B\ K (ii L ) = BM , B^ L ) = Bt(jji L ) , (3.22) 

Bf hh (» L ) = B 2 (» L ), Bl Lh (f, L ) = ^B 2 (ji L ) - 2 -Bz{n L ). (3.23) 

The parameters i?^ LL for all meson systems are also known from most recent lattice simu- 
lations. In this case we use the RG invariant parameters BY LL , usually denoted by Bk and 
£>B q , that are already known from the SM analyses. As these parameters are the same for 
the VRR contributions, we can represent the latter as corrections to the SM box function 
So(xt)- In this manner the VLL and VRR contributions are governed by meson system 
dependent functions 

Si = S (x t ) + ASi (i = K, B d , B s ) (3.24) 

where the ASi are obtained by evolving the Wilson coefficients from ^ down to fif The 
formula for Si is given in (3.32). 

If the unknown 0(a s ) corrections to the Wilson coefficients of non-standard operators 
are assumed to be small, our analysis involves only the values of the coefficients (K), 
P^ R (B d ) and P 2 LR (-B S ) calculated at jXR = 2.5 TeV in the case of box diagrams and hh = 

15 TeV in the case of neutral Higgs contributions. To obtain these values we only need the 
values of B4 and B§. In the case of the K system for /j>l = 2 GeV we have [78] (see table 

16 in that paper) 

£4 = 0.810(41) (31), B 5 = 0.562(39) (46) . (3.25) 

Within the uncertainties we can take the same values for B d and B s systems. In this case 
for fii = fib = 4.6GeV we have [77] 

B 4 = 1.14(3)(6) , B 5 = 1.79(4)(18) . (3.26) 
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with the asymmetric errors corrected using [80] . In the case of box diagram contributions 
(fj, R = 2.5 TeV) we find 

P 2 LR (K) = 73(4) (3) , P 2 LR (B q ) = 4.57(54) (25) . (Box) (3.27) 

For the Higgs contribution (hh = 15 TeV) we find 

Pt R (K) = 88(5)(3), P£ R (B q ) = 5.54(65) (30). (Higgs) (3.28) 

For this calculation we have used the values of m s and nid provided by the lattice-averaging- 
group. We collect them in table 6. The values relevant for VLL and VRR operators are also 
given in table 6 and the final formulae for the mixing observables where all these matrix 
elements enter are presented below. 

3.5 Final expressions for mixing amplitudes 

We now summarise the expressions for the mixing amplitudes Mf 2 , defined in terms of the 
effective Hamiltonian by 

2m Bq (Mf 2 T = (B° q \n A c i= 2 \B° g ) . (3.29) 
We decompose them first as follows 



Then 



where 



M\ 2 = (M? 2 )su + (M? 2 )bh + (A*? 2 )lr = « 2 ) SM + (Mf 2 ) LR . (3.30) 
(M&sm = ^2 F B q B B q m Bq M^ [X^*{B q )] 2 m S* q {B q ) , (3.31) 

S q (B q ) - S ( Xt ) + - [x ^ )]2 . (3.32) 

Here t]b is the known SM QCD correction and t\bI^b ~ 0.95 describes the QCD evolution 
from [iji down to fiw an d is therefore the same for the K° — K° system. So(xt) is given in 
(3.36) below. 

For the LR contribution we first combine the Higgs contributions in (3.15) and (3.18) 

into 

AffiggsCa R (/xr, B q ) = A H oC 2 R (fiH, B q ) + A H +C 2 R (fiH, B q ) . (3.33) 

Then 

c 2 m 2 

(Mf 2 ) LR = -J^fF 2 B m Bq [(A Box C 2 LR (^, B q ))*P£ R fa R ) + (3.34) 
(A Hi gg S C 2 LR ( / iH,i? 9 ))* J P 2 LR ( / UH) 

In these expressions fin = 0{My/ R ) and fin = 0{Mh)- In the case of K° — K° system B q 
should be replaced by K and tjb by r]2- Moreover one should add the known contributions 
from cc and ct box diagrams to (Mj^)sM so that 

(Mf 2 ) SM = ^F 2 K B K m K M 2 w [[X^ L * (K)] 2 m S (x c ) + [A LL * (K)] 2 m S (x t )+ (3.35) 
2Xl L *(K)Xl L *(K) m So(x c ,x t )] , 
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where Fk is the K-meson decay constant and mx the K- meson mass. Here 
Q( \- Q , \ 4x t - llx 2 + x t 3 3x t 3 lnxt 

S {x t ) = S hh {x t , x t ) = — — T2 -^f. , 3.36 

4(l-x t ) z 2(1 - x t y 

S (x c ) = S LL (x c ,x c ) « x c , (3.37) 

3xt 3x 2 In xt 



So(x c ,x t ) = S LL (x t ,x c ) Ri x c 



In 



(3.38) 



x c 4(1 -x t ) 4(1 -xt) 2 

In the last two expressions we have kept only linear terms in x c <C 1, but of course all 
orders in Xf. 

3.6 General anatomy of LR contributions 

In order to see the importance of different NP contributions we rewrite the dominant LR 
contribution in (3.34) by separating contributions of quark mixing matrices from the loop 
integral and QCD running. We obtain 

C 2 M 2 

(Mf 2 ) LR = ^-fFlm Bq ^ij( B q)* R ij( B q) > (3-39) 

i,j=u,c,t 

where we defined 

Rij(B q ) = S LR (x i ,x j ,(3)P^ R (B q ^ R ) 

+ Sl R (x u x 3 ,f3 H )P^ R (B q , hh) (3.40) 
16vr 2 



V2M H G F 



u(s)^Xi(^ H )xj(HH)P2 R {Bq, Hh) 



A ij (B q ) = % R (B q )\f<{B q ). (3.41) 

Here we indicated that the QCD factors P^R are the ones for the B q system. In the 
case of K° — K° system B q should be replaced by K. Choosing M^ R = 2.5 TeV, Mh = 
16TeV • ti(s) 1 / 4 and the central values for the factors Pf given above, we find for the 
matrices R(K) and R(B q ) the results collected in appendix D, for two different choices of 
s. This hierarchical structure of the matrix Rij has an impact on the resulting structure 
of the mixing matrix V R . We discuss this in more detail in section 3.7. 

We observe that for fixed A^ R A^ L the neutral Higgs H° contributions are by far domi- 
nant, followed by gauge boson contributions and rather small charged Higgs contribu- 
tions. This shows that for Mw R > 2.5 TeV the neglect of neutral Higgs contributions for 
masses Mh even as high as 16 TeV, as done often in the literature, totally misrepresents 
the LR story in AF = 2 transitions. In figure 3 we show the relative importance of the 
neutral Higgs contributions to Ru(K) as a function of Mh for different values of Mw R an d 
two different choices for s. We observe that even for M\y R = 400 GeV, which is excluded 
already for many years, and small s the neutral Higgs contributions for Mh < 20 TeV 
account for at least 20% of the total value. For realistic Wr masses M\y R > 2 TeV the 
neutral Higgs contributions become only negligible for masses beyond 100 TeV. For large 
values of s ~ 0.5 the heavy Higgs contribution becomes even more important. As both 
M\y r and Mh are kr dependent decoupling the heavy Higgs contribution would require a 
non-perturbative coupling in the Higgs sector. 
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Figure 3. The relative importance of the neutral Higgs contribution R^ t (see 3.40) to R tt {K) as a 
function of Mjj for different Wr masses for s = 0.1 and s = 0.5 



3.7 Special anatomy of LR contributions 

As can be seen from the formulae in appendix D the matrices R(q) have a very special 
pattern. We first restrict our discussion to the K system. Similar arguments hold for the 
B system but the actual hierarchies change. For now we are only interested in the order 
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10 -6 10 -4 10 - 2 | (3 42) 

10" 4 1(T 2 



of magnitude of the elements and especially in their relative size 

/ KT 9 1(T 6 1(T 4 

R(K) ~ (-1) 

V 

First of all we notice the wide spread of sizes as well as the expected hierarchy towards 
the third generation. The relative shift for changing one generation is about two-to-three 
orders of magnitude. Before summation the elements of R(K) get dressed with elements 
of A(if ) as explicitly given in (3.39). Here we define A (if) as the matrix form of Ay (if). 

Barring fine-tuned cancellations among the various contributions, each of the LR con- 
tributions to AMk and in particular ex has to be suppressed well below the SM contri- 
bution in order to obtain agreement with the data. Thus we can estimate rough upper 
bounds on some of the elements of A (if ) 

|A tt (if)| £ KT 5 , (3.43) 
|A cMc (if)| £ 10- 3 , (3.44) 
|A cc (if)| £ HT 1 , (3.45) 

while no useful bound can be obtained on the remaining elements of A(if ). 

On the other hand for a generic right-handed matrix the matrix A(if ) would exhibit 
the hierarchy implied by the structure of the CKM matrix only 

/kt 1 KT 2 1(T 3N 

1 KT 1 10- 2 I . (3.46) 



|A(if)| 
Recall that 



\1Q- 2 KT 2 10" 



A ij (K) = V£*V%V%*Vk. ( 3 - 47 ) 
Comparing the entries of (3.46) with the bounds obtained in (3.43)-(3.45) we see that 
the ct, tc and tt elements need an additional suppression by one order of magnitude. This 
can only be achieved with the help of an appropriate hierarchy in V^, see (3.47). Explicitly 
we find the constraints 

\V$\\V£\ <; KT 1 , \V$\\VS\ £ HT 1 , \V%\\V*\ & HT 1 . (3.48) 

Following an analogous procedure for the system we obtain 

\V%\\Vi\ ^10- 2 , (3.49) 

while no relevant constraint can be obtained from the B s system. 

Making the plausible assumption that the diagonal elements of V R and in particular 
\V^ \ are 0(1) we obtain 

I^I^HT 2 . (3.50) 

This bound agrees surprisingly well with the one obtained from the explicit numerical 
analysis in section 9.2. We stress however that we have performed here a very rough 
estimate, keeping only the orders of magnitude. It is interesting to see that already this 
naive estimate allows us to understand certain patterns in our numerical analysis. 
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3.8 Comparison of the operator structure in various models 

It is instructive to compare the operator structure in the effective Hamiltonian for 
AF = 2 transitions in the specific LR model considered in the present paper with the 
other specific models containing RH currents considered by us like RSc models [68] or 
some supersymmetric flavour models [81] and in particular in the effective theory approach 
for the RH currents in [44]. To this end let us note that when 0{a s ) corrections at the 
high matching scale are neglected the following dynamics is responsible for the structure 
of the effective AF = 2 Hamiltonian: 

• A tree level exchange of a colourless gauge boson with LH and RH couplings generates 
the operators Q7 LL , Q7 rr and Q\ R - This is an example of Z' models and gauge 
flavour models [82]. 

• A tree level exchange of a gauge boson carrying colour generates the operators Qi LL , 

Qi R an d Qq^- An example is the tree- level exchange of the KK-gluon in RS 



• A tree level exchange of a colourless Higgs scalar generates the operators Q± , Q ^ 
and but as we have seen above at C(e 2 ) only the last operator contributes in 
the model considered. 



• Finally box diagrams with internal charged gauge bosons or carrying both LH 
and RH couplings generate the operators QY LL , Q7 rr and at 0(e 2 ) . 

With this classification in mind it is evident that the effective AF = 2 Hamiltonian at 
the matching scale in [44] corresponds to a tree level exchange of a colourless gauge boson. 
Clearly with RH currents present in this model also the operator QJj ^ s generated through 
box diagrams, but in the presence of QCD corrections this operator is also generated from 
qLR g enera t ec l by the tree level gauge boson exchange in question. An implicit assumption 
in [44] was that this is the dominant mechanism for the generation of Q^ R - 

In the models considered by us there are no flavour changing neutral gauge boson 
exchanges at tree level and the leading mechanism for the generation of are box 
diagrams. The latter can also generate Qi R at 0(e 4 ) but this effect is smaller than the 
QCD mixing generating Q\ R from that we include in our paper. Thus the structures 
of the AF = 2 Hamiltonians considered here and in [44] are in a sense complementary to 
each other. 

3.9 Basic formulae for AF = 2 observables 

We collect here the formulae that we used in our numerical analysis. The mixing amplitude 
M\2 can be decomposed into SM and NP part (i = K, d, s) 



models. 



• A tree level exchange of a Higgs scalar carrying colour generates the operators Q 



iSLL 
1,2 > 





(3.51) 
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and is related for % = K to the relevant effective Hamiltonian through 



2m K (M«)* = (K°\n A e t 2 \K ) , (3.52) 

with analogous expressions for q = d,s. A general formula for the r. h. s. is given in (3.19). 
The Kl — Ks mass difference is then given by 

AM K = 2 [Re (M£) SM + Re (M£) Np ] , (3.53) 

and the CP-violating parameter £k by 

£K = V2(AM K ) cxp [lm (Ml ^ )sM + Im (Mi ^ )np] ' (3 - 54) 

where (p £ = (43.51 ± 0.05)° and k £ = 0.94 ± 0.02 [5, 83] takes into account that (p £ ^ vr/4 
and includes long distance effects in ImTi2 and ImMi2- The value of k £ given here has been 
calculated within the SM using the data on e' je that could also contain NP contributions. 
As analysed in [84] these effects do not have a significant impact on our analysis. For the 
mass differences in the B^ — systems we have 

AM q = 2|(Af? 2 ) SM + (Mf 2 ) Np \ (q = d,s). (3.55) 

Let us then write [85] 

M\ 2 = (Mf 2 ) SM + (Mf 2 ) NP = (Mf 2 ) SM C Bq e 2i ^ , (3.56) 

where 

(M 1 d 2 ) SM = |(M 1 t 2 )s M |e 2 ^ , « 22° , (3.57) 

(Mf 2 ) SM = \(M( 2 )su\e 2i ^ , P s ^ -1° • (3.58) 

Here the phases j3 and (3 S are defined through 

Vtd = \Vtd\e- ip and V ts = -\V ts \e~ 1 ^ . (3.59) 



We find then 



and 



AM q = (AM q ) SM C Bq , (3.60) 



= sin^ + 2</> B J , (3.61) 
^ = 8^(21/3,1-2^), (3.62) 

with the latter two observables being the coefficients of sin(AM^t) and sin(AM s t) in the 
time dependent asymmetries in B^ — > ipK$ and B® — > i/xj), respectively. Thus in the 
presence of non- vanishing <p s d and <p b 3 these two asymmetries do not measure (5 and (3 S but 
(j3 + (fB d ) and (\/3 s \ — <pb s ), respectively. At this stage a few comments on the assumptions 
leading to expressions (3.61) and (3.62) are in order. These simple formulae follow only 
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if there are no weak phases in the decay amplitudes for B® — > t/jKs and B® — > tpcj) as is 
the case in the SM and also in the LHT model, where due to T-parity there are no new 
contributions to decay amplitudes at tree level so that these amplitudes are dominated by 
SM contributions [67]. In the model discussed in the present paper new contributions to 
decay amplitudes with non-vanishing weak phases are present at tree level. However, as 
we demonstrate in section 5 these contribution can be totally neglected when calculating 
S^Ks and S^. 

Now in models like the LHT model and SM4, the only operators contributing to 
the amplitudes and Mf 2 are the SM ones, that is QY LL [67, 69]. Consequently the 
new phases ipB d and ifB a have purely perturbative character related to the fundamental 
dynamics at short distance scales. The situation in the LR model in question is different. 
As now new operators contribute to the Mf 2 amplitudes, the parameters Cs q and ipB q in 
(3.56) are complicated functions of fundamental short distance parameters of the model 
and of the non-perturbative parameters Bi present in Pf"{K) and P^B^). Thus the test 
of the LR models considered with the help of particle-antiparticle mixing and related CP- 
violation is less theoretically clean than in the case of new physics scenarios in which only 
the operator Q^ LL contributes. 

Next, we give the expressions for the width differences AT q and the semileptonic CP- 
asymmetries Ag L 



Ar g / AM q 



r 9 V 



exp 



Re(i¥f M ^ + Im^^ SMsin2 ^ 



MZJ C Bq \M q 12 J C Bq 

r? 2 \ SM cos2^n ^ /r? 2 \ SM sin2v9 B , 



(3.63) 



A ^Hm -^r-HwJ ^f 1 - (3 - 64) 

Theoretical predictions of both AT q and A^ L require the non-perturbative calculation of 
the off-diagonal matrix element T\ 2 , the absorptive part of the B® — B® amplitude as well 
as perturbative QCD calculations. The latter are known at the NLO level [86-90]. The 
most recent results read [10, 91] 

/ F d \ SM / ps \SM 

Re VMV = -5.3(10) -10- 3 , Re h^f) = -5.0(10) -10- 3 , (3.65) 
Im MM = -4.1(6) -10- 4 , Im -if = 1.9(3) • 10" 5 . (3.66) 



\M*J ' ' \M( 2/ 

Finally, we recall the existence of a correlation between and that has been 
pointed out in [97] and which has been investigated model-independently in [98] and in 
the context of the LHT model in [67]. This correlation follows analytically from (3.62) and 
(3.64) when the SM phase (3 S and the first term in (3.64) are neglected: 

/ r q \SMq 

A, ^H^i) k- (3 - 67) 

In [91] it has been pointed out recently that the accuracy of a similar correlation that uses 
AAd s and Ar s [99] instead of Tf 2 and Mf 2 is very poor both for small and large NP phase 
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observable 


experimental value 


SM prediction 


<Ps = -2(P 8 + <p B .) 

Ar s 


€ [-1.04,-0.04] (CDF [92]) 

-0.55±g;|« (DO [93]) 
+0.13(18)(7) (LHCb [94]) 
0.011(37) [95] 
0.075(35)(l)ps- 1 (CDF [92]) 

0.1631°;°!! ps- 1 (DO [93]) 
0.123(29)(8)ps- 1 (LHCb [94]) 


-0.0363(17) [3] 

0.0042(8) [10] 
0.087(21) ps" 1 [10] 


Ad 


-0.12(52)% [96] 
-1.8(11)% [96] 
-0.79(20)% [96] 


-0.041(6)% [91] 
0.0019(3)% [91] 
-0.020(3)% [91] 



Table 1. Theoretical and experimental values of a number of observables related to B s & — B a t i 
mixing. 

(fB 3 - The approximate formula (3.67) is instead very accurate for large S^. In order to 
improve the accuracy also for small values of this asymmetry, in our numerical analysis 
as in our previous analyses in the context of other extensions of the SM, we find such 
correlation numerically by using (3.62) and (3.64) without making any approximations. 

3.10 Summary 

In summary, in this section, we have calculated the NP contributions in the LR model 
in question to the amplitudes and M\ 2 - We have then given formulae for AMk, 
AAig, ek, S^k s , Sij,fa AT q and A§ L in a form suitable for the study of the size of the NP 
contributions. The numerical analysis of these observables is presented in section 10. While 
particle-antiparticle mixing in LR models has already been discussed in the literature, 
our analysis goes beyond these papers as in addition to the full renormalisation group 
analysis and inclusion of all important effects, we search for correlations between various 
observables that have not been studied by other authors. Most importantly, our philosophy 
in performing phenomenology differs from the one used in most papers. Instead for looking 
for bounds on the Wr and Higgs masses we investigate whether the LRM can solve certain 
anomalies present in the flavour data while being consistent with electroweak precision 
tests and the data for tree level charged currents. Moreover, we search for the oases in the 
large space of parameters in which the matrix V R takes special forms that are dictated by 
the data. 
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4 The decays B — > X S)d 7 



4.1 Preliminaries 

The B —7- Xg'j decay in a model with SU (2)^ x SU(2)r x ?7(1) symmetry has been analysed 
by many authors in the past [54-62]. There are basically two classes of contributions: 

• First, the ones resulting from the mixing between Wl and Wr that imply RH cou- 
plings of the SM to quarks. In the SM the LH structure of these couplings 
requires the chirality flip, necessary for b — > sj transition to occur, only through the 
mass of the initial or the final state quark. Consequently the amplitude is propor- 
tional to rrifj or m s . In contrast in LR models the RH couplings allow the chirality 
flip on the internal top quark line resulting in an enhancement factor mt/nn, of the 
NP contribution relative to the SM one at the level of the amplitude. This is the 
contribution mostly studied in the literature. The relevant LO QCD corrections have 
been analysed for the first time within the effective field theory framework by Cho 
and Misiak [57] and have been checked since then by many authors, in particular by 
Bobeth et al. [100], where also NLO QCD corrections to the matching conditions 
at fin have been calculated. In what follows we adopt their results but include NP 
corrections at the LO, while taking into account the known NNLO corrections within 
the SM. 

• The second contribution comes from charged Higgs exchanges. Although in the LR 
models the masses of are O(kr) and numerically significantly larger than few 
TeV, as pointed out in [58] and also analysed in [59-61], the corresponding amplitude 
is also enhanced by nit/nib in contrast to the MSSM where it is proportional to mj 
or m s . Moreover, it does not suffer from the suppression through Wl — Wr mixing 
as is the case of the gauge contribution. As we will see for charged Higgs masses even 
above 10 TeV this contribution cannot be neglected and in fact it can be dominant 
for certain ranges of parameters. This should be contrasted with AF = 2 processes 
where it is generally subleading. 

In the next two sections, we summarise the results for the Wilson coefficients of the 
dipole operators for these two classes of contributions at the relevant matching scales for 
the SM and NP. Subsequently we include RG QCD corrections to these coefficients and 
present the final formula for the branching ratio for the B —> X s j decay. We also present 
the formulae for the CP-averaged branching ratio of the B — > decay and direct CP- 
asymmetries in both decays. 

4.2 Gauge boson contributions 

Adopting the overall normalisation of the SM effective Hamiltonian we have 

AGv 

n cS (b sj) = --^-V t * s V tb [C 7l (fi b )Q 7l + C 8G {fi b )QsG] , (4.1) 

where yu& = 0(mb)- The dipole operators are defined as 

= v^mbSaa^PRbaF^ , Q 8G = ^m h s a a^ P R T^G% . (4.2) 



- 23 - 



In writing (4.1) we have dropped the primed operators that are obtained from (4.2) by 
replacing Pr by Pl- In the SM the primed operators (RL) are suppressed by m s /mb 
relative to the ones in (4.1). As we demonstrate below they can also be neglected in the 
LRM discussed by us. 

The coefficients Cj(/^) are calculated from their initial values at high energy scales 
by means of renormalisation group methods. Before entering the discussion of QCD cor- 
rections we describe here our treatment of LR contributions at high energy scales. We 
first decompose the Wilson coefficients at the scale = (D(M\y) as the sum of the SM 
contribution and the NP contributions: 

Ci(jM W ) = Cf M ( m ) + A LR Q( m ) (4.3) 

and similarly for the primed coefficients. For the SM coefficients we have 



3x? - 2xj , -8xf - bxi + 7x t 

hixjH —, 

2A(x t - l) 3 

-rSMy \ _ ~3%t i . ~ x t + 5%t + %X t 



(W) - i{xt _ JjJ tax, + 2i{xt _ 1)3 = C 7l (z ( ) , (4.4) 



C ^>= 4(^1)4 '■".+ 8(„- l)3 - < 4 ' 5 ' 

The expressions for A LR Ci(/iv^) hi the LRM have been found by Cho and Misiak [57]. 
The by far dominant contribution comes from the induced right-handed part of the Wl 
vertex. At the time of the work of these authors the expected values for M\y R were of the 
order of several hundred GeV and the scale in the LR contributions could be chosen to 
be nw With llr 3> fiw one has to take the effect of large logarithms log(/XR/ ' Liw) into 
account. While a complete RG analysis would be more involved, in the present paper we 
take such effects only approximately into account by simply declaring the result in [57] to 
be valid not at fj,yy but lir. A more involved analysis will be presented elsewhere. 

Adapting the formulae of Cho and Misiak [57] to our notations we find then for the 
LR contributions 



A LR C 77 (/^) = A tb 
A lr C 8G (vr) = A tb 



3x? - 2x t , -5jc? + 31x t - 20 

inx t + 



2(1 -x t ) 3 12(1 -x t y 



-3xt —x? — xt — 4 



lnx t + 



2(1-^)3 1 4(1 -x t f 



(4.6) 
(4.7) 



where 



A tb = m M 2 e m[S] +0(e 4 ). (1.8) 

The Wilson coefficients of the primed operators can be obtained from (4.6), (4.7) by re- 
placing A tb with (A ts )*, where 

A ts = sce 2 e «* (%)+ 0(e 4 ) . (4.9) 
m b \Vt L sJ 

We observe that they are also enhanced by mt/mb in contrast to the primed operators in 
the SM. We stress that mt(mt) and m&(/^&) should be used here. We now give arguments 
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that in order to obtain the leading 0(e 2 ) corrections to the branching ratio we only have 
to keep the contributions in (4.6) and (4.7) while neglecting the contributions from primed 
operators and the LL and RR contributions from NP. Here the LL stands for pure Pl cou- 
plings in the weak gauge boson-quark couplings and analogously for RR with Pi replaced 
by Pr. The unprimed LR contributions result from the Pr coupling in the vertex contain- 
ing the 6-quark, while the corresponding primed contributions result from the Pr coupling 
in the vertex containing the s-quark. This is evident from the couplings A tb and A ts given 
in (4.8) and (4.9), respectively. Now, the LR contributions presented above have been 
obtained by including only SM internal Wl and top quark exchanges taking into account 
the right-handed couplings of Wl in the vertex containing the 6-quark. These contribu- 
tions are important due to the factor mt/rrn, that is absent in LL and RR contributions. 
Moreover, this LR contribution is 0(e 2 ) at the amplitude level and interfering with the SM 
contribution gives also 0(e 2 ) contribution to the branching ratio for B — > X s j. Concerning 
the primed LR contribution to the decay amplitude while being of the same order in e as 
the unprimed LR coefficients and also enhanced by mt/mb, it does not interfere with the 
SM contributions and consequently enters the branching ratio at the e 4 level. Therefore 
it should be neglected for the sake of consistency. Only in the case of a special hierarchy 
of the elements of matrix could this suppression be compensated by the last factor in 
A ts . However, then also other 0(e ) contributions to the rate would have to be included, 
which is beyond the scope of this paper. Finally, as already stressed by Cho and Misiak the 
diagrams with internal Wr exchanges give negligible contributions and a similar remark 
applies to LL and RR contributions from NP as one can easily check. Note that in this 
case the LL and RR contributions are governed by SM loop functions in (4.4) and (4.5) 
which are strongly suppressed when xt = ml/My^ is replaced by xt = m ^/M^ R . 

4.3 Charged Higgs contributions 

In presenting the results for charged Higgs contributions we follow [58, 61] adjusting their 
formulae to our notations and overall normalisations and keeping the phase a. In the 
numerical analysis we set a = 0. As in the case of the gauge boson contributions one 
can demonstrate that the primed operators can be neglected. We also neglect the 
contribution to Cgc (//#). The dominant contribution at bin is given as follows [58] 4 



A H± C 77 (/i„) = -u(s) \sc^ a A 1 H+ (y)+2s 2 c 2 A 2 H+ (y) 



'VI 

where the function u{s) has been defined in (3.17) and [58] 



3y 2 - 2y 5y 2 - 3y 

3(1-0)3 n 2/ + 6(1 



(4.10) 



(4.11) 



A 2 H+ (y) = -A SM (y) - A l H+ (y) , (4.12) 



A SM (y)= -2Cf 7 %), 2/=^jr- (4-13) 



4 In [58] tan/3 = k/k' — c/s has been used. 
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The last function is given on the r.h.s of (4.4). In these contributions we set m% = mt(//ff). 
For large s close to the limit s — > 1/V2 a strong enhancement of the contribution 
through the factor u(s) is possible. As we discussed in section 2.4 taking this limit, cor- 
responding to k' = k, is phenomenologically not viable. Interestingly, the inspection of 
gauge boson and charged Higgs contribution shows that provided the element V^ R has only 
a small phase, these contributions always enhance the branching ratio for B — > X s ~f which 
brings the theoretical value in (4.23) closer to the data in (4.24). 

4.4 QCD corrections 

In order to complete the analysis of B — > Xg'j we have to include QCD corrections which 
play a very important role in this decay. In the SM these corrections are known at the 
NNLO level. In the LR model a complete LO analysis has been done by Cho and Misiak 
[57]. On the other hand Bobeth et al provide matching conditions to the Wilson coeffi- 
cients of LR operators at hr relevant for a NLO analysis, that is 0(a s ) corrections to the 
coefficients (4.6) and (4.7). However, performing a complete NLO analysis would require 
additional complicated calculations and in view of many parameters present in this model, 
such an involved analysis is certainly premature. In view of these remarks we proceed as 
follows: 

• For the SM contribution we use the full result at the NNLO [101]. 

• For the LR contribution we use explicit LO formulae which we obtained on the basis 
of [57] and the recent paper [102]. We set the fib scale to the one used for the SM 
contributions, this means fib = 2.5 GeV. While the NNLO SM contribution is not 
sensitive to this choice, some sensitivity is present in the LO LR contribution but in 
view of several parameters involved this uncertainty is not essential. We are aware 
of the approximate nature of this treatment of QCD corrections in the NP part but 
we think that such an approach is sufficient before the discovery of Wr. 

Thus the basic formula for C 77 (//?,) used by us reads: 

C 77 (/x 6 ) = C 77 (^) SM + A LR C 77 (/x 6 ) + A H± C 7y (fx b ) , (4.14) 

where 

A LR C 77 (/i fe ) = k 7 (jmr) A LR C 77 (/i fi ) + ks(jmr) A LR C SG (iu R ) + A cb KLR (fu R ) , (4.15) 
A H± C 7l (fi b ) = k t (hh) A h± C 77 (^) . (4.16) 
Here A LR C 77 (/x K ), A LR C 8G (w?) and A H± C 77 ( / u ff ) can be found in (4.6), (4.7) and (4.10), 



respectively. A cb is given as follows 



m b V cb 



A cb = p sce 2 e ia ^. (4.17) 



The term proportional to A cb is absent in (4.16) as, being related to the mixing with new 
charged current operators, it is already included in (4.15) and above the scale /ir this 
mixing does not take place. 



- 26 - 





1 TeV 


2.5 TeV 


10 TeV 


15 TeV 


K7 


0.457 


0.427 


0.390 


0.380 




0.125 


0.128 


0.130 


0.130 


KT T! 


0.665 


0.778 


0.953 


1.005 


ro 


0.504 


0.475 


0.439 


0.429 


01 R 


-0.052 


-0.043 


-0.025 


-0.019 


k 7 


0.857 


0.801 


0.731 


0.712 


k 8 


0.044 


0.060 


0.078 


0.082 


KLR 


0.063 


0.099 


0.156 


0.175 


P~8 


0.874 


0.824 


0.760 


0.743 


PLR 


-0.033 


-0.044 


-0.056 


-0.058 



Table 2. The NP magic numbers for AC??- and ACg<9 at fi b = 2.5 GeV and fi t (m t ). 



Finally, k's are the NP magic numbers listed in Tab. 2 [102], calculated taking a s (Mz = 
91.1876 GeV) = 0.118 5 . 
For later purposes we also give 

A LR C 8G (p b ) = P8 (pr) A LR C 8G (/x i? ) + A cb PLR (^ R ) , (4.18) 

with the NP magic numbers pi listed in Tab. 2. 

4.5 The branching ratio 

For the branching ratio we follow the strategy of [102] which used the results of [103]. One 
has then 

Bt(B -> X s7 ) = Br(B -> X sl )su + ABr, (4.19) 



where 
with 
Next 
and 



ABr = R 



2Re(C 77 (^) SM A LR C 77 ( Atfe )) + |A LR C 77 ( Atfe )| 2 J , (4.20) 
A^CVyG^) = A LR C 77 (/i 6 ) + A H± C 77 (^) . (4.21) 
R = 0.00247 , C 7y (p b ) SM = -0.353 , (4.22) 



Br(B -> X s7 ) S m = (3.15 ± 0.23) x 10" 4 (4.23) 

are extracted from [103]. We refer to [102] for details. In particular in obtaining the 
value for C 77 (^) SM we have taken the non-perturbative corrections into account. Strictly 
speaking the last term in (4.20) is of 0(e 4 ) and should be dropped together with the 
contributions of primed operators. We recall that experimentally [95] 

Bt(B X^)^ = (3.55 ± 0.26) x 10" 4 . (4.24) 

so that ABr > is favoured implying A LR C 77 (fib) < 0. This is guaranteed for ReV t R > 0. 



3 We thank Emmanuel Stamou for providing this table. 
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4.6 Dissecting the LR contributions to Br(B — >■ X s -f) 

Let us next calculate the fraction of contributions of the full NP contribution at the 
level of the branching ratio. As seen in (4.20) and (4.21), if we neglect the last term in 
(4.21) this is simply given by A H± CV 7 (/ift)/A LR CV 7 (//&). In figure 4 we plot this ratio 
as a function of Mh for different values of M\y R and two choices of s. For these plots 
we set y R to be the identity matrix, making the formulae V R independent. Since the 
main contributions from charged Higgs and gauge bosons are both proportional to this 
assumption has a very limited impact. 

As in our previous discussion of AF = 2 contributions in section 3.6, the Higgs ex- 
changes (this time H^ 1 ) again turn out to be very important. For realistic Wr masses 
above 2 TeV and H± masses below 20 TeV, the contributions account for at least 20% 
(s = 0.1) or 50% (s = 0.5) of the total NP effect. As expected for the large value of s = 0.5 
the contribution is dominant. 

4.7 The decay B -> X d j 

Of considerable interest is also the decay B — > for which the measured CP-averaged 
branching ratio reads [95] 

(Bi(B -> Xtf)) = (1.41 ± 0.49) x 10~ 5 , (4.25) 

to be compared with the most recent SM value [36] 

(Br(B -> X dl )) SM = (1.54±g;§f) x 10~ 5 . (4.26) 

While the data agree with the latter estimate there is still significant room for NP con- 
tributions. Yet already these results can imply important constraints on the extensions of 
the SM. Indeed, it was recently pointed out by Crivellin and Mercolli [36] that the contri- 
butions of primed operators to (Br(5 — > Xjpf)) could be significantly constrained implying 
a bound on \V3\. This bound is approximately 3.5 times stronger than what is found for 
the best- fit solution in [44], where this decay has not been considered. In the notation of 
the present paper the bound in [36] reads 6 

sce 2 \V$\ < 1.4 x 10~ 4 . (4.27) 

As our numerical analysis shows the electroweak precision tests imply see 2 ^ 1 • 10~ 3 . 
Therefore as long as \V^\ < 0.14 the bound in [36] is satisfied in our paper. Moreover as 
we can see in section 10 the typical values for \V^\ that are consistent with ex are below 

R 1 0.047 (._06) 
1 td] ~y 0.113 (s = 0.1) v 1 

when no fine-tuning constraint is taken into account. Otherwise the allowed ranges for \ V^\ 
become even smaller. For further details see section 9.2. In all cases \ V^\ is well below the 



Note that this bound does not take into account the charged Higgs contribution. 
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Figure 4. The relative importance of the charged Higgs diagrams in the NP contributions to 
Br(B — > X s j) for s = 0.1 (top panel) and s = 0.5 (bottom panel) as a function of the Higgs mass 
for different Wr masses. 

bound in [36], since e.g. a large s implies small values for \V^\. As we explained before, 
in our model the contributions of primed operators are 0(e 4 ) and consistently neglected in 
our analysis. While a detailed analysis of C(e 4 ) contributions to Br(B — > X^j) would be 
required to assess the numerical importance of primed operators in this case, the discussion 
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a 




< 


a 77 


T 


a 7e 


088 


a 8 


7.8221 


0.4384 


-1.6981 


0.8161 


4.8802 


-0.7827 


0.0197 


0.5680 


a 8e 


a 87 


a\ 


4 


4 


a 87 


4, 




-0.0601 


0.1923 


0.3546 


-0.0987 


2.4997 


-0.0487 


-0.9067 


-0.0661 



Table 3. The relevant a,; parameters from [105]. 



above indicates that in our model they are well below the NP contributions of unprimed 
operators when all additional constraints are taken into account. 

Clearly, the improved data on the CP-averaged branching ratio in (4.25) and the finding 
of [104] that the CP-averaged branching ratio contains only small hadronic uncertainties, 
invites us to consider this decay as well. In analysing it we use, as done in [36], the 
formulae of [105] . In doing this we should emphasise that these formulae strictly speaking 
apply only to cases where the RG evolution below the jxw scale is the same as in the SM. 
Due to the presence of new LR current-current operators, additional contributions to the 
Wilson coefficients of the dipole operators are present. They are represented by the last 
term in (4.15). However, in the model considered by us these contributions are very small 
and can be neglected. We find then 



(Br(B -> X dl )) = ^ 



V td V tb 



2 

[a + P 7 + P 8 + P 78 ] , (4.29) 



where N = 2.57 x 10~ 3 and 

a = a + a ee \ed\ 2 + a r e Re(ed) , (4.30) 

P 7 = a 77 \R 7 \ 2 + a r 7 Re(R 7 ) + a r 7e Re(R 7 e* d ) , (4.31) 

Ps = a 8S \Rs\ 2 + a r 8 Re(R 8 ) + a r 8e Re(R s e d ) , (4.32) 

P 78 = a r 87 Re(R s R* 7 ). (4.33) 
The values of the parameters are collected in table 3 and 

^ 7 -^W)' ^-^w' (434) 

where the C 7 j(fj,t),C 8 G(tH) denotes the total Wilson coefficients including both the SM 
contributions without QCD corrections, as given in (4.4) and (4.5) and the NP contributions 
as in (4.15) and (4.18) but evaluated at /it and not 



A LR C 77 (/i t ) = R 7 (v R ) A LR C 77 (^) + R^r) A LR C 8G (/i J? ) + A cb K LK (fj, R ) , (4.35) 

A H± C 77 (/i 4 ) = k 7 { m ) A H± C 77 (/i H ) , (4.36) 

and 

A LR C 8G (^) = ~ P8 { m ) A LR C 8G (fi R ) + A cb ~ PljK ( m ) . (4.37) 
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The NP magic numbers hi and pi are listed in table 2. Finally 



(4.38) 



For — > X s 7 one should just replace d by s. 
4.8 CP asymmetries in B -> ^s,d7 

A very sensitive observable to NP CP violating effects is represented by the direct CP 
asymmetry in b — > S7, i.e. ^cp(fr - ► $7) [106], in particular as the perturbative contri- 
butions within the SM amount to only +0.5% [105]. The corresponding asymmetry in 
b — > d'j transition is much larger but could in principle also provide a useful test. These 
asymmetries are defined by q = (s, d) 



r(g -> Xfy) - T(B -» X ql ) 
T(B -> X g - 7 ) + T(B -> X ?7 ) 



(4.39) 



and have been studied in particular in [107, 108] and more recently in [105], where further 
references can be found. In the context of the LR models these asymmetries have been 
analysed in [60]. Unfortunately, a recent analysis [109] shows that these asymmetries, 
similar to other direct CP asymmetries suffer from hadronic uncertainties originating here 
in the hadronic component of the photon. These uncertainties lower the predictive power of 
these observables and in the case of b — > S7 the authors conclude that only if experimentally 
Acp(b —> 57) was found below —2%, one could consider it as a signal of NP. In order to 
get a rough idea whether in the models considered by us the perturbative part could be 
affected strongly by NP, we use the formulae in [105], which are compatible with [107, 108]. 
The formulae with hadronic contributions that are rather uncertain can be found in [109]. 
We have then 

2 

lm.{a\R 7 + a\R& + a\e q + a l m R^R* 7 + a\ e R 7 e* + a l 8e R 8 e* q ) 
(Br(B ^ X ql )) 

(4.40) 

with the values of a« collected in table 3. 



Ac?{b ->■ qj) 



100 



T/L* T/L 

V tq V tb 



5 Constraints from tree level decays 
5.1 Preliminaries 

In this section we address the constraints from tree level decays. The possible new tree 
level contributions arise from the new right-handed couplings of the Wl gauge boson, from 
the exchange of the heavy Wr gauge boson, and from the heavy charged Higgs boson. 
For a detailed consideration of leptonic and semi-leptonic decays it is necessary to have a 
closer look at the lepton sector within the LR model. In this context in order to derive the 
Feynman rules we use the findings of [63, 64]. After transformation to mass eigenstates the 
light neutrinos are dominated by their left-handed contribution with a small right-handed 
admixture, while the heavy Majorana neutrinos are mainly given by the right-handed 
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Br(vr nv) = 0.9998770(4) [112] 


/tt = 129.5(17) MeV 

£ l (X\ n ClKQA (A A \ 

J+yJ) = U.yoo4(44J 
f K /f w = 1.1931(53) 
f Ds = 248.9(39) MeV 
f B = 205(12) MeV [114] 


J + {v)\Vus\ — 0.21b3(5J 
f K /U\V us /V ud \ K ^ = 0.2758(5) [113] 


Br(D s -> tu) = 0.0529(28) 
Br(B -»• tu) = 1.64(34) • 10" 4 
F{l)\V cb \ B ^ D * lu = 0.03604(52) 
G(l)\V cb \ B ^ mu = 0.0423(15) [95] 


F(l) = 0.908(17) [115] 


G(l) = 1.074(24) [116] 



Table 4. Values of the most important experimental and theoretical quantities used as input 
parameters for the constraints on tree level charged currents. 

neutrinos again modified by a small left-handed contribution. This mixing angle can be 
constrained by the masses of heavy and light neutrinos. Assuming reasonable masses for 
the light neutrinos and the heavy neutrinos not to be lighter than 100 GeV, the Yukawa 
couplings have to be very small and cause this mixing to be at most O(10~ 6 ) [110]. In 
agreement with [111] we find these mixing effects to be negligible. This is in particular 
a good assumption taking into account that they have to compete with effects of 0(e 2 ) 
being roughly of 0(1O -3 ). Furthermore tree level decays will only take place into the 
light neutrinos. As the Wr coupling to light neutrinos appears first at 0(e 2 ) and the 
Wr propagator yields another 0(e 2 ) suppression factor, we conclude that for the lepton 
couplings to charged gauge bosons at C(e 2 ) only the SM couplings are relevant for the tree- 
level decays in question. Concerning the tree level contributions of charged Higgs bosons 
it is sufficient to include only leading order couplings. In [44] a detailed analysis of the 
constraints on the elements of the matrix (denoted there by V") and implications for 
the CKM matrix V L have been presented. As the charged Higgs contributions to tree level 
decays in the present model turn out to be negligible (see below), basically all the results 
obtained in [44] can be taken over by making the following identification: 

e L = 0(e 4 ), e R = cse 2 , V = e ia V R , (5.1) 

where the quantities on the l.h.s are the ones used in [44]. For completeness we summarise 
these results in our notation, extending the discussion of the impact of the right-handed 
currents in tree level decays on the mixing induced asymmetries S^k s an d and on the 
CKM phase 7 as extracted from B — >■ D decays within the SM. We use our findings from 
[44] with updated numerical values. 

5.2 Constraints on the mixing matrices V 1, and V R 

In what follows we use the notation 

\Vij\v = |V# + cse ia e 2 V$\ , = \Vj - cse ia e 2 V%\ , (5.2) 

for the combinations of left- and right-handed contributions entering vector and axial vector 
couplings. As Gf enters tree level decays we note that the correction to the width of the 
//-decay in this model are at C(e 4 ) and have no impact on our estimates of V R . 
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5.2.1 d 

The study of super-allowed + — > + transitions yields the constraint on the u — > d vector 
current [112] 

\V ud \ v = 0.97425(22) . (5.3) 

A constraint on the axial component can be obtained from the pion decay ir + — > fJ- + v, 
including the known radiative corrections [117], with the result 

\VudU = 0.981(13) . (5.4) 

5.2.2 u ->• s 

The vector current s — > u transition can be obtained from K — > itlv decays [113] with the 
result 

\V us \y = 0.2257(12) . (5.5) 

The axial s — > u transition can be constrained by combining K — )■ [iv and tt — > [iv. The 
result reads 

\V us \a = 0.2268(32) . (5.6) 

5.2.3 c d 

The constraint 

\V cd \ v = 0.229(25) (5.7) 

can be obtained from the decays D — > Klv and D — > ixtv [112]. \V cd \ can also be obtained 
from neutrino and anti-neutrino charm production off valence d quarks. As interference 
terms between left- and right-handed quarks are suppressed by m d , we can safely neglect 
them. We can therefore directly apply the constraint to left-handed c — > d transitions and 
get [112] 

|^|= 0.230(11). (5.8) 

Due to the large uncertainties these data do not provide significant constraints on the LR 
parameter space but we include them in our numerical analysis. The situation should 
change in the future when more precise data and lattice inputs will be available. Similar 
comments apply also to c — > s transitions discussed below. 

5.2.4 c-^s 

From semileptonic D decays and D s — > t + v one finds 

1^ = 0.98(10), \Vc\a = 0.978(31) . (5.9) 

5.2.5 b^u 

This is a place where the RH currents could enter in a potentially important manner. If 
the SM value for \ V u b\ is determined from the inclusive semileptonic mode B — > X u £u, the 
interference term between LH and RH contributions is totally negligible, so that the SM 
result carries over to the LR model. We have [112] 

\V^ b \ = 4.27(38) • 10" 3 , (5.10) 
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which should be regarded as the true value of this CKM element. The vector b — > u 
transition can be probed by B — > n£u. This implies the constraint [112] 

\V uh \ v = 3.38(36) • 1(T 3 . (5.11) 

Finally the axial b — > u coupling is determined from Br(B — > tv) implying the con- 
straint 

\V ub \ A = 4.70(56) -1(T 3 . (5.12) 

5.2.6 6->c 

Similar to the case of b — )■ u transitions, the inclusive determination of \V c b\ in the SM 
carries over to the LR models, where we have 

\V%\ = 41.54(73) • 1(T 3 . (5.13) 

On the other hand from B — > Dlv and B — > D*iv transitions one finds the constraints 

\V ch \ v = 39.4(17) • 1(T 3 , \V cb \ A = 39.70(92) • 10~ 3 . (5.14) 

The good agreement between the two exclusive determinations suggests a small contribu- 
tion from right-handed currents and consequently LR effects cannot explain the tension 
between inclusive and exclusive determinations. 

5.2.7 t^-d,s 

The t — > d and t — > s transitions cannot be measured from tree level decays, so that no 
constraint is obtained on Vj^ ts and V^ ts . 

5.2.8 t b 

In the SM \Vtb\ can be determined from the ratio of branching ratios Br(t — > bW)/15i(t — > 
qW) = \ Vtb\ 2 assuming unitarity of the CKM matrix. As interference terms are suppressed 
by m q /m t (q = d,s,b), this method can directly be applied to LR models. The best 
available constraint stems from DO [118] and reads 

|^|=0.95(2). (5.15) 

This measurement is 2.5a below the SM expectation Vi ~ 1. Since as we see later the tight 
tree level constraints on other elements on the CKM matrix do not allow for a significant 
deviation of from its SM value, this tension persists in the LRM. The improved data 
from the LHC will tell us whether this is a real problem. 

5.3 Charged Higgs contributions 

The charged Higgs couplings are given by combinations of the Yukawa couplings of the 
fermions participating in the interaction. Therefore for all processes not involving the top 
quark the charged Higgs contribution is very small and turns out to be negligible in all 
cases. The only decay where this is not completely obvious is B + — > t + v t due to the chiral 
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suppression of the W + contribution. However it turns out that the H + contribution is 
negligible also in that case. Indeed, the Higgs tree level diagrams appear with a suppression 
factor of m 2 B / M'jj and hence in the LR models in which is very heavy this contribution 
is negligibly small in comparison to the effects from right-handed charged gauge boson 
couplings discussed above. 

5.4 A comment on Z — >■ bb 

Even if the electroweak precision tests are considered in the next section it is useful to 
discuss already here the impact of new contributions on the Z —> ff couplings. Denoting 
the effective diagonal couplings of the Z to down-type quarks as follows 

£fff = — (9t + g« $ R <f<? R ) Z, , (5.16) 

we find (see appendix) at tree level 

= (-\ + \s 2 w) ~ l~/ R c 2 R e\ (5.17) 

{9%) = \s w - C f{cl + \sl)e\ (5.18) 

where the first terms are the SM contributions. The experimental determination of the 
effective couplings of the Z bosons to light down quarks resulting from the global fit of 
electroweak data collected by the LEP and the SLD experiments is in a very good agreement 
with the SM. For b quark couplings one finds [119] 

(g b L )exp = -0.4182 ± 0.0015 , (5.19) 
(sgjexp = +0.0962 ± 0.0063 . (5.20) 

While the result for the LH coupling is consistent with the SM prediction, there is a large 
disagreement between data and SM expectation in the RH sector: 

(A<$) exp = G$) cxp - Q$) SM = (1.9 ± 0.6) x 10~ 2 , (5.21) 

It is evident from the formulae above that with e 2 = 0(1O -3 ) as required by our analysis 
of electroweak precision constraints the corrections to Z — > qq are marginal and that while 
not spoiling the agreement of the SM with the data for the remaining couplings, these 
corrections cannot help in removing the anomaly in question. In this context we would like 
to emphasise that the flavour independence of the corrections would preclude the solution 
anyway without spoiling the agreement for the remaining couplings. Finally, the sign of 
corrections being strictly negative would make the disagreement with the data even worse. 

5.5 Impact on CP asymmetries 

5.5.1 CP asymmetry in B — > DK and the angle 7 

In the SM the angle 

7 = &rg(-V ud V: b /V cd V; b ) , (5.22) 
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can be measured by means of CP asymmetries associated with tree level B — > DK decays. 
The theoretically cleanest method has been proposed in [120]: Studying B^ — > DK^ 
together with a Dalitz plot analysis of a following multibody D decay allows for determining 
the weak phase 7 without prior assumptions on the hadronic parameters involved. Belle 
[121] obtained 

7 = 77(11)°, (5.23) 

where we combined the several uncertainties in quadrature. The details of the Dalitz 
plot analysis can be found in [120] - as they are relevant mostly for the extraction of the 
hadronic parameters from the data we do not repeat them here. The measurement relies 
on the interference between the two modes 

A(B- D°K-) = A B , (5.24) 
A(B- -> D°K-) = A B r B e i( - SB -^ a) , (5.25) 

with relative absolute value r B , relative strong phase 5 B and relative weak phase 7 e g. In the 
SM the two amplitudes are governed by the operators {cb)y _ a(su)v - a and (ub)v-A(sc)v-A 
mediated by the W boson, respectively. Therefore in this approximation 7 e fj = 7. 
In LR models there are two new contributions to B^ — > DK^- decays: 

1. The small RH coupling of W leads to contributions of {V — A) ® (V + A) structure. 
Unless the corresponding hadronic matrix elements are strongly enhanced with re- 
spect to the standard {V — A) ® (V — A) ones, this contribution is subleading with 
respect to the one obtained from W exchange as it is suppressed by an additional 
factor sc < 1/2. We neglect this contribution in view of the presently large uncer- 
tainties in the experimental determination of 7, keeping in mind that it should be 
reconsidered once precise data become available. 

2. Similarly to the case of S^k s (see below) we also have the contributions of Wr with 
purely right-handed couplings, that generate the (V + A) Cg) (V + ^4) operators. As 
QCD is non-chiral the relevant matrix elements are equal to the SM ones. The LR 
contribution considered here is therefore purely perturbative. 



The interfering amplitudes then get modified as follows 



m 2 w v*(y*f 



A(B- -+ D K~) = A(B- -> D"K-) SU 1 + r " , (5.26) 

V Wr V cb V us J 

1V1 W R V ub V cs J 

Here we neglected the QCD running of the LR contribution from Myy R to M^/. Then the 
relative weak phase is given by 



7eff = 7 + arg 



\ + r M w fv£(y*y v*(v*r 



Wr V v cb v us v ub v cs 



(5.28) 
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5.5.2 The impact on S^Ks and 

In the SM the decay B d — > J/ipKs is dominated by one single decay amplitude arising 
from a tree level W exchange and giving rise to the operator (bc)v-A(cs)v~A- Using the 
conventions of [122] the weak phase of the decay amplitude A{B d — > J/ipKs) is to a very 
good approximation given by 

<£ D = axgV^Vi;~0. (5.29) 
Together with the phase of — B^ mixing 

M = argV^~-/3 (5.30) 

the time dependent CP asymmetry therefore measures 7 

S^ Ka = sin[2(0 D - <j> M )} = - sin[2 arg = sin 2/3 , (5.31) 

cb cs 

In LR models we have new contributions to the b —> ccs tree level decay. First of all the 
W now has a small right-handed coupling proportional to e 2 ~ O(10 -3 ). Second, also Wr 
exchanges contribute at the e 2 ~ 0(1O~ 3 ) level. Assuming that the QCD matrix elements 
do not significantly affect the hierarchy of new contributions, the direct Wr exchange is 
dominant, in particular if s <C 1. As QCD is a non-chiral theory, we find for the matrix 
elements governing the decay in question 

(J/iliKs\(bc) V -A(c8)v-A\B d ) = (J/i>Ks\(bc) v+A (cs) v+A \B d ) . (5.32) 

Hence the LR correction to the decay amplitude can be expressed in terms of short-distance 
physics only. We find 

/ M 2 v r* v r\ 

A(B d -> J/^K S ) = A(B d -> J/V^s)sm 1 + xr J f » . (5.33) 

\ 1V1 W R V cb V cs ) 

Here x encodes the RG running from Myy R down to M\y and r = s^/c^Sr = 9rI 9l- 
Consequently 

S^k s = sin[2(/3 + p Bd + 5(j) D )] , (5.34) 



with 



/ M 2 V R*yR\ 

iD ---A i+xr it«m)' (5 - 35) 



We estimated x ~ 1.1. 

The same arguments apply to B s — > J/ip(f), hence we find 

S^4 = - sin[2(/? s + ip Bs + Hd)] ■ (5.36) 

Note that 8<p£) is universal in both transitions. 



7 For the estimates of the uncertainties in this relation see [123]. 
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5.5.3 Summary 

Using the expressions above we find that the impact of new contributions on the determi- 
nation of 7 and of new phases ipB 3 and cpB d is negligible. This can be traced back primarily 
to Wr being much heavier than Wl- 

6 Electroweak precision constraints 

6.1 Preliminaries 

Electroweak precision tests (EWPT) provide very strong constraints on basically any model 
of NP. The SM agrees to a high accuracy with the measured values of roughly 40 low-and 
high energy observables [112, 124, 125]. Very often such tests are performed with the help 
of the famous S-T variables [126], which test the allowed size of the oblique corrections, 
but in order to get the full picture including the LEP II data the study of all observables 
should be favoured. 

Since a full analysis of the EWP data is clearly beyond the scope of this work we 
use the results of [35] where a full analysis of EWP observables in a number of models 
with SU(2)i x SU (2)2 x U(\)x gauge symmetry has been performed. The model denoted 
there by LR-T is precisely the model studied in our paper. Most useful for our study are 
tables IX and X, where the authors of [35] provide formulae for corrections to the most 
constraining observables in a given model in terms of the fit parameters. These formulae 
allow then to find allowed regions in the space of the fit parameters that are consistent 
with the data on EWP data. This analysis is independent of the flavour parameters of the 
matrix V^. On the other hand the analysis of flavour physics observables is affected by 
the choice of EW parameters. This is in particular the case for s, to which the 6—7-57 
decay and to a certain extent also AF = 2 processes are very sensitive. Because of its 
enhancement by charged Higgs contributions 6—7-57 requires s < 0.64. In this sense we 
are able to separate the EWP analysis from the flavour analysis. We first determine the 
allowed ranges for EW parameters and then perform the flavour analysis for a certain EW 
benchmark point, varying only the parameter s. 

In what follows, we very briefly describe the main points of [35], simultaneously ad- 
justing their notation to ours. 

6.2 Basic structure of the analysis 

As seen in (2.23) and (2.24) except for g s , nit, Mh and the phase a, which we set to 
zero without loss of generality, we have six model parameters that are relevant for this 
section. For the present analysis it is convenient to make a change of variables and to fix 
the following input {reference) parameters 



which are most precisely measured. In what follows we use, as in [35], the value of the 
QED coupling constant in the MS scheme: 



G F 



M z 



a e = a{M 2 z ) 



(6.1) 



l/a(M§) = 127.916 ±0.015. 



(6.2) 
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In the case of the fit parameters that parametrise NP contributions, we choose 

e, s R , s, (6.3) 
which are related to the ones in [35] through 

x = , c^ = cr, sin 2/3 = 2sc . (6.4) 

The goal of [35] is to find the allowed ranges for the three parameters in (6.4) by fitting 
37 EWP observables to the existing data. To this end the SM contributions are included 
at the tree and one-loop level, while NP contributions are taken into account at the tree 
level. This is sufficient for our purposes. Before describing the manner in which we use 
the results of [35], let us make a few comments on how our choice of reference parameters 
affects the 0(e 2 ) corrections in our formulae that have been exclusively written in terms of 
model parameters. The latter parameters have been distinguished in [35] through a tilde. 
Expressing these parameters through the parameters in (6.1) and (6.4) we find that 

• Gf and M%, as expected, do not receive any 0(e 2 ) corrections. This can be explicitly 
verified by expressing the model parameters in terms of the parameters (6.1) and (6.4) 
and inserting in the formula for Mz given in appendix B. The 0(e 2 ) correction cancels 
when v and sw are defined by 8 

v 2 - 1 s 2 c 2 - ™ {M2z) (6 5) 

v - ^mc~ F ' SwCw ~ v2m 2 g f • (6 - 5) 

• On the other hand Mjy, M z > and Myy R are given as follows: 

M w = M z c w 



,2 2 4 



1 + 2 c 2 - s 2 { 4 ' 
w w 



(6.6) 



c W b R c R c W b R 
where we have dropped 0(e 2 ) corrections in the last equation, s\y and cyy are defined 
through (6.5) and the one-loop SM corrections have not been shown in (6.6). 

The implications for our analysis are as follows: 

• The formula (6.6) after the inclusion of SM electroweak loop corrections together 
with other observables considered in [35] can be used to constrain the fit parameters 
in (6.4). However in the present paper we simplify this analysis by requiring that 
the shift in Mjy due to NP contributions is consistent within 2a with the difference 
between its measured value and its SM prediction 

(AM w f p = e 2 ^^v4 " 2sV ) = M w ~ M ^ ( 6 - 8 ) 
= 0.036(34) GeV, 
where the SM and measured value are taken from [124, 127]. 



3 Our v is by y/2 smaller than the one used in [35] 
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• On the other hand (6.7) can be used, as done in [35], to find lower bounds on the 
masses of these new gauge bosons. In turn incorporating the corresponding bounds 
from collider experiments can put additional constraints on (6.4). 

• Finally, in case and Z' will be discovered and their masses precisely measured, 
the formulae in (6.7) will allow to take these values as reference parameters reducing 
the number of fit parameters to one. 

6.3 Basic constraints 

Using the dictionary in (6.4) and tables IX and X of [35] we find the following expressions 
for the four observables of interest: 



fhad the partial branching fraction of Z — > qq with the NP correction 

<5cThad/0"had,SM = 



2 4 

-1.13^ - 0.142^ + 0.0432(2s 2 c 2 ) 
4 4 ; 



the forward-backward asymmetry ^4fb(^) 



SA FB (b)/A FBtSM (b) 



c 2 c 4 
-30.0^ + 67.6^ 
4 4 



2„2\ 



20.6(2s^c 



the weak charge Qw(Cs) of the caesium-133 nucleus 



5Qw(Cs)/Q w ,sm(Cs) 



-0.855^ - 0.145(2s 2 c 2 ) 
4 ' 



the left-handed coupling measured in deep inelastic i/-N scattering 

S(g^) 2 /(9L,SU? = [0-0219 + 0.478c| + 0.210c 4 - 1.42(4s 2 c 2 )] e 2 



(6.9) 



(6.10) 



(6.11) 



(6.12) 



A more detailed description of the observables considered here is presented in e.g. [35, 128]. 
The correct treatment of the constraints is then given by the following formula 



|EXP - SM(1 + CON) | < V(AEXP) 2 + (ASM(1 + CON)) 2 , 



(6.13) 



where EXP and SM stand for the experimental and SM value, respectively. CON stands 
for one of the conditions we have listed above. Furthermore we use the constraint from 
the W boson mass as given in (6.8), where the SM and experimental values are taken from 
[124, 127]. 

We also incorporate the direct experimental constraints from muon decay of the TWIST 
Collaboration [129]. They present two bounds on the ratio of gn/gi times the mixing an- 
gle and one on gL/dR with respect to the Wr mass. As s\y/(cwSr) = gn/gi holds, we 
can make further simplifications as in the first constraint the factor gn/gi cancels. In our 
notation their results read then 



see 1 < 0.020 



and 



cyySR 
sw 



M Wr > 578 GeV 



(6.14) 
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at 90% C.L.. 

Moreover we take into account the direct experimental bounds on Mw R . While [130, 
131] report stringent constraints > 2 TeV on the mass of heavy W gauge bosons, these 
bounds only apply to the case of very light right-handed neutrinos which escape detection. 
In case of a non-negligible M VR < Mw R the constraints on the masses have been analysed 
in [132], with the bound on M\y R reaching up to 1.7 TeV depending on the right-handed 
neutrino mass, with only 240 pb -1 of data analysed so far. On the other hand if M VR > 
M\y r the branching ratio of Wr decaying to leptons is e 2 suppressed and the dominant 
decay mode is into quarks. In that case the searches for resonances in the dijet mass 
distribution [133, 134] apply, yielding the constraint M\y R > 1.5 TeV. We note that all 
these bounds depend on the ratio ghldR through the production cross-section for Wr. In 
order to keep our analysis independent of the details of the LRM neutrino sector, we will 
assume Mr > 2 TeV throughout our analysis. 

Finally for completeness we include the bounds g R < Att and g' 2 < Att in order to 
guarantee perturbativity of the gauge couplings. 

7 Strategy for the numerical analysis 

7.1 Preliminaries 

Having at hand all the relevant formulae for AF = 2 processes, B — > Xd )S l observables, 
EWP observables and tree level decays in the LRM in question, we are ready to perform 
a numerical analysis taking all existing experimental constraints into account. 

The first question which one could ask is whether there are regions in the parameter 
space of the LRM for which all observables can be found in their allowed ranges. As 
the LRM has a large number of parameters, in particular in the mixing matrix y R , it is 
probable that there are such regions even for an LHC accessible mass M\y R . 

The next interesting question is if there are experimentally allowed simplified parametri- 
sations for the RH mixing matrix with only a few free parameters. There are a few popular 
yet already excluded matrices of this kind. We however derive inspiration from those and 
propose a new and simple matrix. 

In general regions in the parameters space fulfilling 

• the fine-tuning of parameters is small, 

• the anomalies in the present data are softened or removed, 

are particularly interesting. These regions can be considered as different 'oases' in the 
parameter space which are disconnected and lead in general to different phenomenology. 

7.2 Anomalies in the flavour data 

Before proceeding with outlining our strategy we now review the known flavour anomalies. 
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7.2.1 The ek — S^Ks anomaly 

It has been pointed out in [5, 84] that the SM prediction for ek implied by the measured 
value of S^k s = sin 2/3, the ratio AM^/AM, and the value of \V c b\ turns out to be too 
small to agree well with experiment. This tension between ek and S^k s has been pointed 
out from a different perspective in [4, 8, 135, 136]. These findings have been confirmed by a 
UTfit analysis [2]. The CKMfitter group having a different treatment of uncertainties finds 
less significant effects in ek [9]. Indeed taking the experimental value of S^k s = 0.679 ± 
0.020, \V c b\ = 0.0406, the most recent value of the relevant non-perturbative parameter 
Bk = 0.737 ±0.020 [137-141] resulting from unquenched lattice calculations and including 
long distance (LD) effects in ImFi2 and ImMi2 in the K° — K° mixing [5, 83] as well as 
recently calculated NNLO QCD corrections to ek [142, 143] one finds [143] 

|ex|sM = (1-81 ± 0.28) -10- 3 , (7.1) 

visibly below the experimental value l^/flexp = (2.228 ± 0.011) • 10" 3 . On the other hand 
sin 2/3 Rj 0.85 ±0.05 from SM fits of the Unitarity Triangle (UT) is significantly larger than 
the experimental value S^k s = 0.679 ± 0.020. This discrepancy is to some extent caused 
by the desire to fit e k [4, 5, 8, 84, 135, 136] and Br(£> + ->■ t + v t ) [8]. For the most recent 
discussions including up to date numerics see [10, 11, 144]. As demonstrated in [5, 84], 
whether the NP is required in ek or S^k s depends on the values of 7, \V u b\ and \ V c b\. The 
phase 7 should be measured precisely by LHCb in the coming years while \V u b\ and \ V c b\ 
should be precisely determined by Belle II [145] and SuperB [146-148] provided also the 
hadronic uncertainties will be under better control. Here we concentrate briefly on \V u b\- 

7.2.2 I Kb I problem 

As already mentioned previously there is the tension between inclusive and exclusive de- 
terminations of \ V u b\ with the exclusive ones in the ballpark of 3.4 • 10 -3 and the inclusive 
ones typically above 4.0 • 10~ 3 . As discussed in [41] an interesting possible solution to this 
problem is the presence of RH charged currents, which selects the inclusive value as the true 
value (see our discussion in section 5), implying again sin 2/3 ~ 0.80 [44]. Unfortunately in 
the LRM the increased value of the Wr mass precludes this explanation if only points with 
acceptable fine-tuning are considered. Indeed the factor see < 10~ 3 from EWPT data 
and the bound on \V^\ from FCNC processes as found by us, see section 10.3 for details, 
imply that the effect of RH currents is simply too small. Therefore one cannot state that 
the inclusive value is the favoured one from the point of view of the LRM and both the 
exclusive and inclusive determinations are equally valid. 

7.2.3 Possible solutions 

As discussed in [4, 5] and subsequent papers of these authors a negative NP phase ipB d 
in B® — B^ mixing would solve the ek — S^k s anomaly, provided such a phase is allowed 
by other constraints. Indeed, this is evident from (3.61). With a negative <fiB d the true 
sin 2/3 is larger than S^k s , implying a higher value on \ek\i in reasonable agreement with 
data and a better UT-fit. This solution would favour the inclusive value of \V u b\- On the 
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other hand as stressed in [5] a sizeable constructive NP physics contribution to ek would 
not require an increased value of sin 2/3 relative to the experimental value of S^k s an d the 
exclusive value of \V u b\ would be favoured in this case. Clearly also in such an analysis 
AM^s should be considered. 

7.2.4 Enhanced value of S^^ 

This topic became rather hot in the last years but the situation is still unclear at present. 
The Tevatron data for combined with the results for the same sign dimuon asymmetry 
of DO indicated last spring non-standard CP violation in the -B^-system corresponding to 
~ 0.8 [10-12]. On the other hand the subsequent measurements of CDF and DO 
[92, 93] and the first more accurate results for from LHCb [94] indicate consistency 
with the SM. Yet, the range 

-0.1<S^<0.4 (7.2) 

summarised recently in [149] leaves still a lot of room for sizeable NP contributions. Let us 
hope that the future data from Tevatron and in particular from the LHCb, will measure 
this asymmetry with sufficient precision so that we will know to which extent NP is at 
work here. 

7.3 Addressing the flavour anomalies in the LRM 

Now we want to discuss how the LRM can in principle address the anomalies described 
above. As the LRM has a large number of parameters, in particular in the mixing matrix 
F R , it is obvious that for certain choices of these parameters, most, or even all of these 
anomalies can be removed. Yet, as we will see the situation is more involved than one might 
think at first sight. In order to better understand this matter we introduce an analytic 
expression for 1/ R valid in a specific oasis of the parameter space. This matrix is given in 
terms of only few parameters. Other oases will be considered later on. 

In this manner we see that already the present data give us some hints on possible 
structures of the V R matrix. The fact that the low energy flavour data can give us such 
information about physics taking place at scales of few TeV is clearly remarkable and com- 
plementary to direct collider searches for NP at the LHC, where it is basically impossible 
to learn about the structure of V R . With improved flavour data we should be able to find 
out which of the oases considered by us, if any, is favoured by nature. 

7.3.1 Various scenarios for \V u b\ 

In view of the fact that the LRM cannot solve the \ V u b\ problem, as mentioned above and 
discussed in section 9.2, we introduce two scenarios for \V u b\- In particular, we investigate 
the correlations between ek, S^k s , S^, AM^ and AM S in these two scenarios for \V u b\ 
setting 7 = 68°. These scenarios are defined as follows: 

1. Small | Kb | 

In this scenario we set 

\V ub \ =3.4x 10- 3 . (7.3) 
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In this scenario within the SM we find S^k s ~ 0.675 in agreement with the data but 
£k ~ 1.8 x 10~ 3 is visibly below the data. 

2. Large \V ub \ 

In this scenario we set 

\V ub \ =4.4x 1(T 3 . (7.4) 

This is the case, considering again the SM, in which ek is consistent with the data 
while S^k s ~ 0.82, significantly above the data. 

We now illustrate with an example how all these discrepancies can be solved with a 
special form of the matrix that depends only on four parameters: S13, S23 and two 
phases (f>i and ^2, all chosen to be in the first quadrant. We want to stress that these 
scenarios are only used in our discussion regarding the simplified matrix introduced in the 
following section. As already discussed in section 5 in our general analysis we treat the 
constraints on \V u b\ and 7 like any other measurement. 

7.3.2 Understanding anomalies through V R 

In this section we want to propose a new simplified RH mixing matrix and subsequently 
illustrate how it can be used to solve the flavour anomalies. This matrix has been found 
with the goal to have very simple expressions for the AF = 2 observables while satisfying 
all existing constraints. To this end it turned out useful to set S12 = in as in this 
manner, as seen in (2.27), five entries in this matrix simplified considerably. Moreover, the 
AMk constraint has been relaxed in this way. The reduction of the number of phases to 
two turned out also to improve the transparency. 

The resulting special form of V R is then given by 



-else"* s 13 \ 

C23Sl3e~ t(t>1 -S23e"^ 2 C23C13 / 



(7.5) 



In order to derive simple formulae we restrict the mixing angles to the phenomenolog- 
ically viable ranges 

5i3^0.02, 523^0.2, (7.6) 

with the constraints coming from the and B s system, respectively. Without further 
assumptions on the phases involved, a much stronger combined constraint 

S13S23 £ 10~ 5 (7.7) 

can be derived from ek- In this range of parameters we find then to a very good approx- 
imation that the terms proportional to A^ R A^- L dominate as far as AM^, ek, S^Ks and 
are concerned. For AMk also A^ R A RL is relevant but the appearance of the same 
phase 4>\ in and makes its contribution to ek to be negligible. We find then 

G 2 M 2 

(Im M£) LR = \R tt (K)\ x -J^F 2 K m K \Vt d \\Vt\c23hzhz sin(0 2 - fa - P + ft) . (7.8) 
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For — B® mixing we have 



C Ba e 



Using the matrix (7.5) we obtain 



\Rtt{B q ) 



S (x t )B Bq r]B 



(7.9) 



K H (B d ) 
.^ L (B d ). 

and consequently 



\V L \ 
I v td\ 



R 



X? R (Bs) 
^ L (B S ) 



\vh\ 



(7.10) 



sin 2ip Bd 
sin 2(p Bs 



\Zd\ C23S13^ 



\Zs\ S 23 V£ 



sin(^>i — j3) 



Here 



Cb s \V ts] 
Rtt{B q ) 



\V L 
I v td, 

yR 

Yj~ sin(0 2 - (3 S ) 



(7.11) 
(7.12) 

(7.13) 

S (xt)B Bq r]B 

Setting then Myy B = 2.5 TeV, s = 0.1 and M H o = M H ± = 16TeV we find then for central 
values of other parameters 



1^(^)1 = 9.1 



\Rtt(B q 



0.57. 



\ z d\ 



0.36, 



0.34. 



(7.14) 



These quantities summarise the dominant Slr contribution that governs the box contri- 
butions with (Wl,Wr) and (Wl,H^) as well as neutral Higgs tree level exchanges. As 
seen in (3.40) it includes also QCD corrections, hadronic matrix elements and depends on 
the masses of Wr and heavy Higgs particles. If only the box contributions involving Wl 
and Wr were taken into account the NP contributions to ek would be reduced univer- 
sally by 85% allowing for larger values of the mixing parameters and phases in these three 
expressions. The Higgs dominance can be read off figure 3. 

The fact that three quantities in question depend only on two new phases fa and fa 
implies interesting correlations in the two scenarios for \ V u b\ that we illustrate now. These 
correlations depend of course on the two positive real values of S13 and S23 as we see below. 

• In scenario 1, S^k s within the SM agrees well with the data implying that fa ~ (5. 
This in turn implies that ek and are governed by the phases fa — 2/3 and fa, 
respectively. The desire to reproduce the experimental value of ek requires fa > 
2(3 — f} 3 , and can be enhanced simultaneously. 

• In scenario 2, ek in the SM agrees well with the data implying that fa ~ fa — (3 + f3 s . 
Consequently S$k s and are this time governed by the phases fa — 2/3 and fa, 
respectively. The desire to reproduce the experimental value of S^k s requires again 
fa > 2/3 — (3 S which also leads to an enhanced value of S^. 
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Figure 5. Correlation between and Ek in scenario 1, for fixed 4>\ = and §13 = 3 • f 3 . The 
central SM value is indicated by the red dot. 



In summary the phases 4>\ and 4>2 should satisfy 




P, 



(scenario 1) 
(scenario 2) 



(7-15) 



> 2/3 - p s . 



(7.16) 



As we can see there are relations between the NP phases in the and B s systems but 



In figure 5 we show the correlation between ek and in scenario 1. To this end we 
set 4>i = (3 and S13 = 3 • 1CT 3 and vary §23 < 0.2 and 2/3 — (3 S < < ft + 2/3 — /3 S , imposing 
the experimental constraints on AM S and AM^. The SM point is reached for §23 = 0. The 
black curves correspond to constant 4>2 values as indicated. We observe that for not too 
small S23 ^0.05 and 45° < 4>2 < 50° good agreement with the data for ek can be achieved 
while at the same time enhancing w. r. t. its tiny SM value. For larger values of (f>2 
a smaller S23 is required in order to fit the ek data. At the same time the effects in 
become much smaller. Within scenario 1 this simple structure for V R can thus solve the 
ek anomaly. The measured value of S^^ can then be used to measure the values of S23 
and <j>2- 

In figure 6 we show the correlation between AM^ and AM,, normalised to their SM 
values, after imposing the constraint from ek- Again <j>\ = f3 and §13 = 3 • 10~ 3 are fixed, 
while S23 < 0.2 and 2/3 — f3 s < 4>2 < vr + 2/3 — f} s are varied. Since the observed values 
for both mass differences lie below their SM predictions, a suppression is welcome also in 
this case. We observe that AM^ is very close to its experimental central value and almost 
constant, since it depends to very good approximation only on 4>\ and §13 which are fixed. 
Also AM S can easily be suppressed and brought closer to the data, which in turn favours 
moderate enhancements of S^^ ^ 0.4 (light-blue points in the figure) or very large effects 
S^tj, ~ 1 (red points) - the latter being disfavoured by the LHCb data. 

As a specific example for how the anomalies are solved in scenario 1, we quote the 



they are more involved than the ones characteristic for 2HDM- 



MFV 



[75, 76]. 



L.l 

2 1.0 

S 

< 

S" 0.9 



0.8 
0.7 

<: 

Figure 6. Correlation between AMj and AM, in scenario 1, for fixed 0i = f3 and §13 = 3 • 10 -3 , 
and imposing the constraint from ek- 

following parameter point: 

5i3 = 3 • 1(T 3 , s 23 = 0.03 , 0i = p = 21° , 2 = 50° . (7.17) 
With this choice of parameters we find 

\e K \ ~ 2.2 • 10" 3 , 5^ « 0.27, (7.18) 

and 

AM d ps 0.51 ps" 1 , AM S ps 17.4, ps" 1 (7.19) 

barring in mind the theoretical and parametric uncertainties that have been omitted here 
for the sake of simplicity. 

Let us stress that this solution to the flavour anomalies is in fact the simplest one 
possible within scenario 1. Setting any of the four parameters in (7.5) to zero would spoil 
the relations discussed above. 

Let us now turn our attention to scenario 2, in which ek is in good agreement with 
the data, but S^k s needs to be suppressed w. r. t. its SM value « 0.82. We show now that 
although this simple structure of can yield some improvement over the SM situation, 
the tensions cannot be fully resolved in this case. 

To this end we show in the left panel of figure 7 AM^/ (AM^)sm as a function of S^k s 
for different values of Sis < 0.01 and 2/3 — j3 s < 02 < 150° and fixed 0i = — ft + /3 S . We 
observe that in order to simultaneously fit the data on S^k s an d AM^ 02 is required to 
lie roughly in the range 90° . . . 120°. 

In the right panel of figure 7 we show AM S / (AM s )sm as a function of for different 
values of S23 < 0.1 and 2(3 — f3 s < 02 < 150°. Here we observe that in order to obtain 
a suppression of AM S w. r. t. the SM prediction, as favoured by the data, a much smaller 
phase 02 ^ 70° is required. An enhancement of over its SM value is then obtained 
automatically. Such a low value for 02 would in turn allow only for a mild suppression of 
SipKs with respect to its SM value. This exercise shows very clearly how important it is 
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Figure 7. left: AAf ( i/(AAf ( j)sM as a function of S^k s f° r different values of S13 < 0.01 and 
2/3 — fi s < (j)2 < 150°. right: AM s /(AM s )sm as a function of for different values of §23 < 0.1 
and 2/3 — (3 S < <fi 2 < 150°. The SM values are displayed as black dots, while the experimental la 
regions are indicated by grey bands. 

to consider various observables simultaneously. Each of the panels in figure 7 taken alone 
allows to remove the anomalies in and B s systems, respectively. But when they are 
taken together a clear tension between the suppression of AM S and suppression of S^Ks is 
present. 

In view of the present non-perturbative uncertainties, it is not possible to rule out 
this structure of V R as a solution to the flavour anomalies within scenario 2, however we 
note that this scenario is disfavoured since only a slight amelioration of the SM tensions is 
possible. On the other hand, as seen before, this simple matrix works beautifully in solving 
all anomalies within scenario 1. 

In fact the conclusion is much more general: Within scenario 2 it is not possible to 
solve the flavour anomalies if the tt contribution dominates all observables in question. To 
see this, let us parametrise 

Vu = \Vi\e- i ' t> \ V£ = -\v£\e-^ , V$ = \Vi\e~^ . (7.20) 

The constraint from ek in (7.8) again implies <f>\ ~ 02 — /3 + fi s - In the formulae for Bd tS 
mixing ((7.11), (7.12)) $2 is then replaced by 4>2 — 4>3- However, since the same combination 
of phases appears in both meson systems it cannot help to ameliorate the tension between 
Bd amd B s data. Thus in order to solve the flavour anomalies in the LRM in case of a 
large \V u b\ value, charm contributions have to be relevant. 

These correlations are reminiscent of the ones found in the context of a 2HDM with 
flavour blind phases (2HDM MFV ) [75, 76]. The difference is in the direct contribution to 
ek m scenario 1 which was absent in the 2HDM in question. 

We should also emphasise that with being close to unity, in both scenarios the 
branching ratio Br(S — > X s ^f) is also enhanced over the SM value bringing the central 
theoretical value of this branching ratio closer to experiment. 
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While showing this example we do not claim that the matrix in (7.5) is the only one 
that is capable of removing simultaneously (scenario 1) or ameliorating (scenario 2) all 
anomalies in a correlated manner, but possibly this is the simplest matrix achieving this 
goal. One can also check that moving to different quadrants of the mixing angles and 
phases would generally imply different correlations between the observables in question 
and this fact illustrates how one can determine V R once the data on S^, S^,k s and \ V U \\ 
improve. 

7.4 Nominal input parameters 

We have collected the input parameters required for the numerical analysis in tables 5 
and 6. The values of other parameters like Pf and the tree level constraints on the CKM 
matrix have been given in previous sections. At this stage it is important to recall first 
the theoretical uncertainties in the constraints used by us. Considerable progress in lattice 
calculations has been made in recent years reducing the uncertainties in Fb s and Fg d and 

also in yj Bb s Fb s and J Bb a Fb a down to 5%. This implies an uncertainty of 10% in AM^ 
and AM S within the SM. The numerical values and their respective errors are given in 
table 5. Even bigger progress has been made in the case of the CP-violating parameter ek, 
where the decay constant Fk is known with 1% accuracy. Moreover the parameter Bk 
is known within 3% accuracy from lattice calculations with dynamical fermions [137-141] 
and an improved estimate of long distance contributions to ek reduced this uncertainty 
down to 2% [5, 83]. Also the calculation of NNLO QCD perturbative corrections to ek 
[142, 143] improved the theoretical status of ek by much. The situation with other Bi 
parameters describing the hadronic matrix elements of other operators is much worse as 
clearly seen in the errors of the factors P" in (3.27) and (3.28). Here a significant progress 
is desired. 

Finally let us recall that the CP-asymmetries and S^k s have rather small hadronic 
uncertainties 9 and the hadronic uncertainties in the ratio AMd/ AM S amount to roughly 
3%. Also the theoretical uncertainties in the rate B — > X s ^ decay are below 10%. 

On the experimental side AMd, AM S and ek are very precisely measured, so that their 
experimental errors can be neglected for all practical purposes, while S^k s is known with 
an uncertainty of ±3%. AMjc, while very accurately measured, is subject to poorly known 
long distance contributions and we only require that (AMx) exp is reproduced within ±30%. 
Finally the rate for the B — > X s ^y decay is known within the accuracy of 10%, comparable 
to the theoretical one. 

8 A numerical (re-)analysis of EWPT constraints 

In this section we revisit the constraints from EWPT following the analysis in [35] and 
adding a few new constraints as outlined in section 6. 

9 See [123, 150] for more details and references therein. 
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= 1.16637(1) • 10" 5 GeV" 2 
M w = 80.399(23) GeV 
a(M z ) = 1/127.9 
a s (M z ) = 0.1184(7) 
sin 0w = 0.23116(13) 
mP K = 497.614(24) MeV 
AM K = 0.5292(9) • lO^ps" 1 
|ejr| = 2.228(11) • 10~ 3 
m Bd = 5279.5(3) MeV 
m Ba = 5366.3(6) MeV [112] 


7/1 = 1.87(76) [143] 


773 = 0.496(47) [142, 151] 


7/2 = 0.5765(65) [152] 


rjB = 0.55(1) [152, 153] 


Fk = 156.0(11) MeV 
B K = 0.737(20) 
F Bd = 205(12) MeV 
F Bg = 250(12) MeV 
B Bd = 1.26(11) 
B Bs = 1.33(6) 


AM d = 0.507(4) ps- 1 
AM S = 17.77(12) ps- 1 

Tb s — i-Ai J-i^Oj PS 

T Bi = 1.519(7) ps 

sin(2/3)^ c5s = 0.679(20) [95] 


F Bd ^B Bd = 233(14) MeV 
F Bs ^B Ba = 288(15) MeV 

r>B s /£>B d - i-.VO{<) 

£ = 1.237(32) [114] 


m c {m c ) = 1.268(9) GeV[114, 154] 
m t (m t ) = 163(1) GeV 
m 6 (2.5GeV) =4.60(3) GeV 





Table 5. Values of the experimental and theoretical quantities used as input parameters. 





2GeV 


4.6GeV 


172GeV 


2.5TeV 


15TeV 


m n (/j)(MeV) 


2.09(0)(9) 


1.74(6) (7) 


1.15(8)(5) 


0.97(8)(4) 


0.88(8)(4) 


m d (/i)(MeV) 


4.73(0)(11) 


3.94(1)(9) 


2.61(2)(6) 


2.19(2)(5) 


2.00(2)(5) 


m s (/i)(MeV) 


93.6(2)(11) 


77.9(3)(9) 


51.6(4)(6) 


43.4(4)(5) 


39.5(4)(5) 


m c (/j)(MeV) 


1089(7)(0) 


907(6) (0) 


601(5)(0) 


505(4)(0) 


460(4) (0) 


m b (fx)(GeV) 




4.074(19) (0) 


2.702(14)(0) 


2.268(12)(0) 


2.068(12)(0) 


m t (/x)(GeV) 






162.3(10)(0) 


136.3(9)(0) 


124.2(8)(0) 



Table 6. The NLO running quark masses at different scales. The first and the second parenthesis 
shows the statistical and systematic error, respectively. 



8.1 Numerical procedure for the EWP tests 

As discussed in section 6 we can reduce the free parameters in the electroweak sector down 
to three. We randomly scan over the full ranges of < s < l/y/2, 0.1 ^ sr < 1 and 
< e < 0.1. For every parameter point we evaluate the x 2 function corresponding to the 
constraints discussed in section 6. 

8.2 Electroweak precision constraints 

First of all we want to show the correlations between the three electroweak parameters 
we have chosen in section 6, that is e, s B and s and in this manner establish the allowed 
ranges for these parameters to be used in the analysis of FCNC processes. In figure 8 
we show the correlation of s and e while encoding ln(x 2 /nd. .f.) as a colour spectrum. As 
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one can see smaller values of e are clearly favoured by the data as indicated by the colour 
gradient. Additionally one can see an area of a slightly darker blue at e ~ 0.03 and non- 
trivial s. This area allows to soften the disagreement with Ap B as measured at LEP. This 
measurement has to be taken with a grain of salt because the competing SLD experiments 
did not measure a departure from the SM. We also note that values of e above 0.08 are 
disfavoured. 

In figure 9 we show sr as a function of e. A preference for smaller values of e can again 
be observed. For values of e ~ 0.03 the value of sr has to be non trivial and even above 0.6 
in order to keep x 2 low. As in figure 8 there is a slightly darker region for e = 0.01 . . . 0.03 
and sr = 0.7 . . . 1 where the Ap B constraint would be fulfilled within 2a. 

Figure 10 shows Mw R as a function of e. As expected the mass of the heavy gauge 
boson Wr increases with a decreasing e. The wide spread is due to other parametric 
dependencies. Interestingly in this picture the dark blue region implies an upper limit on 
M\y R . We conclude that in case Ap B as measured at LEP was true we would be able to 
put this model under pressure if My/ R is not found below roughly 7TeV. 

For our analysis of AF = 2 observables we fix two of the electroweak parameters to a 
point with intermediate values 

Si? = 0.80, e = 0.03. (8.1) 

The parameter s, which turns out to be very relevant for flavour observables, is varied in 
the range 

0.1<s<0.6. (8.2) 
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ln(x 2 /n dof ) 
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Figure 9. In this figure we show sr as a function of e. The colour spectrum corresponds to 



ln(x 2 /n d . . f 
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Figure 10. In this figure we show Mw R as a function of e. The colour spectrum corresponds to 
ln(x 2 /™d.o.f.)- 

This range of parameter points is contained in the dark blue area in figures 8-10. In this 
way we are able to incorporate the EWP constraints into our flavour analysis in a quite 
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general manner. Changing the electroweak reference point (8.1) would affect mostly the 
masses My/ R and Mh- The Higgs mass Mjj additionally displays a strong dependence 
on the parameter s. We already studied the effect of changing the masses on AF = 2 
observables and Hr(B —> X s ^f) in sections 3.6 and 4.6, respectively. In terms of more 
tangible parameters our choice corresponds to the masses 

16 

M Wr ps 2.6 TeV , M H rj - TeV , (8.3) 

R y/1 - 2s 2 

where we choose 03 = 8 in the Higgs potential. We emphasise that increasing the value 
Mh by much can only be done at the price of loosing perturbativity in the scalar sector. 
The parameters relevant for our flavour analysis then read 

1.3 < k/k' < 9.9, r=^fw0.48, 8.9 • 10~ 5 < see 2 < 4.3 • 10 -4 . (8.4) 

91 

Note that using the constraints outlined in section 6 we are not able to fulfil Ap B at 
below 1.8(7. This can in part be accounted for by the approximate nature of the constraints 
used but in our opinion this should not be the dominant effect. This suggests that the LRM 
cannot explain the LEP anomaly, but only soften it. 

9 A general study of V R 

In this section we study the structure of the RH matrix V R . In order to perform our 
general analysis of and of flavour observables we varied all 13 parameters in the V L 
and y R matrices in their allowed ranges. For the parameters of this means . . . ir/2 for 
the mixing angles and . . . 2tt for the phases. In the case of V L parameters the situation is 
simpler as due to the smallness of ek and the tree-level constraints one is able to restrict 
the allowed ranges beforehand. We impose all available constraints from AF = 2, Br(S — > 
Xgj) as well as the tree- level measurements of CKM parameters within 2<7. Our choice of 
electroweak parameters is described in section 8.2. 

9.1 Allowed ranges for V R 

First we study the matrix V R in full generality, aiming to identify the valid regions of 
parameter space. To this end we allow the parameter s to vary in the range 0.1 . . .0.6. 
Due to the unitarity of the matrix y R it is sufficient to consider the allowed ranges for the 
absolute values of V^., and in order to obtain a complete picture of the possible 
size of all elements. Therefore in figure 11 we show the correlations between pairs of these 
elements of V^. The value of s is encoded in colour, where all plots show points with 
large s in the front layer. Generally the allowed regions grow with decreasing s. First we 
observe that large regions of the parameter space are already excluded. In addition all three 
elements can reach their maximal values ~ 1 individually, however not simultaneously. The 
two branches of points are easier to disentangle in the 3D version of these plots, as can be 
seen in figure 12 for s = 0.1. 
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Figure 11. |V^|, \VR\ and |V^*| as functions of each other with s encoded in colour 



In order to keep the presentation as transparent as possible, in the following we restrict 
ourselves to showing results only for fixed values of s when necessary. In the course of 
our analysis of AF = 2 constraints we concentrate on the case s = 0.1. For the study 
Br(£? —7- X s ^) large values of s are interesting since they lead to enhanced effects. In a 
few cases we consider the s dependence explicitly. 

9.2 A fine-tuning study 

We now turn to investigating the fine-tuning necessary for a point to fulfil all the experi- 
mental constraints. 

It is well known that models which predict sizeable contributions to the LR AF = 2 
operators can lead to dangerously large fine-tuning of several observables. For example the 
RS model with custodial protection has to struggle with huge effects in ek and a potentially 
large fine-tuning [68, 155]. 

First let us define the measure of fine-tuning 

-. Nob s Wobs / ns^ \ 

Nobs *H N ohs f-' j V °i °Pj J 

1=1 1=1 J 
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Figure 12. The absolute values of \V^ S \, \V^\ and \V^\ with log Agg d encoded as the colour. 



used throughout our analysis. Here Abg is the wei l known Barbieri-Giudice (BG) measure 
of fine-tuning [156]. Note that the sum should only contain observables which are in fact 
fine-tuned in some part of the parameter space. The overall fine-tuning defined in equation 
(9.1) is compatible with the more sophisticated fine-tuning measure proposed by Athron 
and Miller [157]. Since the numerical application of the Athron-Miller fine-tuning measure 
requires much more computing power than the BG measure while not providing additional 
insights we decided to stick with AgQ d . We note that by definition Abg is sensitive only to 
fine-tuning in terms of cancellations between various contributions, but not to accidentally 
small parameters. 

In figure 12 we show the absolute values of |V^|, \V^\ and \V^\ as a three dimensional 
plot for s = 0.1. The colour corresponds to logAgQ 1 and we plotted points with low fine- 
tuning in front of ones with large fine-tuning. This however does not preclude points with 
high fine-tuning to lie in regions dominated by points with low fine-tuning. First of all we 
observe that the constraints allow for three clearly distinct scenarios. 

• In the first scenario \V^\, \V^\ and |V^| are small. The fine-tuning can be very 
small and increases slowly with |V^;|. We call this scenario the "normal hierarchy" 
scenario. 

• We find the second scenario in the case of small and large \ V^\. This scenario 
is not restricted to one corner but allows points along the |T^| axis as well. The 
fine-tuning is in most cases high and increases further with increasing \V^\. Taking 
only points with low fine-tuning this scenario corresponds to the so-called "inverted 
hierarchy" scenario. 
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contributions of the heavy Higgses. 



• For small |V^|, large \V*\ and \V$\ < 0.30 we find the third scenario. This scenario 
exhibits large fine-tuning for all points, hence it is completely eliminated if we require 

A BG d < 10 - We P° int out that while this 

scenario is very fine-tuned it is also not 

rigorously excluded. 

Second we find that while significant regions of the parameter space suffer from large 
fine-tuning, there exist ranges of parameters in which the fine-tuning is small and all 
experimental constraints, in particular the ones from ex and AMk, can be satisfied. 

In the following we restrict ourselves to points in parameter space which exhibit only 
a small level of fine-tuning AgQ d < 10. In this case we find the following limits: 

1^1 < 1.2 .10-' ^ |K5|<{^ (0.2) 

Note that the constraint on |V^| is much more stringent than the one on \V^\. The 
"normal" and "inverted hierarchy" scenarios introduced above are then defined by 

IK* I < 0.3, (normal hierarchy) (9 3) 

\V*l \ > 0.9. (inverted hierarchy) 

The first scenario leads to a hierarchical structure of V R with small off-diagonal elements, 
while the second one inverts the hierarchy of the 2, 3 submatrix giving a very different 
structure. 

Let us now consider how these results change if the heavy Higgs contributions are 
omitted. In figure 13 we show the allowed parameter points obtained from a scan analogous 
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Figure 14. arg(#i(Bg)) as a function of \,^t{B q )\ with on the left-hand side and B 8 on the 
right-hand side. The colour coding corresponds to S^,k s an d respectively. 



to the one that lead to figure 12, taking into account only the gauge boson contributions. 
We observe that now the points cover a significantly larger region of parameter space, even 
if only points with low fine-tuning are considered. This drastic change reflects once more 
the importance of the heavy Higgs particles. 

9.3 A closer look at the "normal hierarchy" scenario 

Now we study in more detail the "normal hierarchy" scenario. Recall that we consider only 
parameter points with small fine-tuning AgQ d < 10. In order to investigate the possible 
contributions of the matrix V R to flavour processes it is useful to study the relative size 
of the dominant LR contribution compared to the SM contribution. From the hierarchical 
structure of V R in this scenario together with the pattern of the matrix R, see appendix 
D, we expect the tt contribution to dominate. Therefore it is convenient to define (q = 
K, B d ,B s ) 

\l R (q)\? L (q) _ A^ R (g) 
\\\q)\\\q) A t LL (g) 
In order to better understand the results obtained in this scenario it is useful to in- 
troduce a Wolfenstein-like parametrisation of the right-handed matrix by expanding in 
A = \V*\. Defining s\2 = A, S13 = BX and §23 = AX and then expanding in A we arrive at 



„g i (0c6-</>tb+0ti;+0iid-<Piis)^ e i i.<t>cb-<t>tb+<t>ts) ^ _ 1 ^2 _|_ j\ ^2^ yle^^A 

^ e**"* ^Ae^'-^ - Be^*"-*""'! A 2 -Ae ilpts \ e**" (l - ^A 2 

(9.5) 

Note that this matrix cannot be reduced to the matrix in V R in (7.5) and is valid in a 
different oasis in the parameter space. In figure 14 we show the correlation of arg (^t(B q )) 
and \^t{Bq)\. On the left hand side of figure 14 we show the situation in the B^ system. 
In order to gain more information from this correlation we encode S^k s as the colour of 
a point. First we notice that the absolute value \^t{Bd)\ only allows enhancements up 
to 1.5. The phase arg(j^£(Z^)) on the other hand is allowed to be in the whole [— tt, tt] 
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range. However there are restrictions on the phase as well. For the maximal allowed 
enhancement of the absolute value the phase is restricted to be close to or ±7r and for 
a phase of ±7r/2 the absolute value is restricted to be below roughly 0.8 and 0.5 for a 
phase of n (—it) respectively. The observed shape reminds us of the model-independent 
constraint on (M^ 2 )np- Since the measured values of AM^ and S^k s are m good agreement 
with the SM prediction (although somewhat on the low side), there is not much room left 
for NP. Furthermore as AM^ suffers from larger theoretical uncertainties than S^k s the 
possible size of the NP contribution is largest if its relative phase is close to or ±7r. 
If the tt contribution was the only contribution to (M{^ 2 )np, its phase would be equal to 
avg(^ t (Bd)) and we would observe strong peaks for these phases in the plot. The small 
"washout" indicates however that other LR contributions cannot be neglected in this case. 

On the right hand side of figure 14 we show the situation in the B s system. In this 
plot the colour corresponds to S^. We immediately see that the situation is completely 
different from the Bj system. The enhancement in the absolute value \^t(B s )\ is allowed 
to be up to nearly 7 while the phase &rg(^t(B s )), though allowed to be somewhere in the 
whole range, is very strongly correlated with the absolute value and S^. For \S^\ close 
to zero the absolute value is restricted to be smaller than 1.5 or bigger than 4.7 while 
the phase is either roughly free or very close to zero respectively. For big \S^\ the phase 
is restricted to two very specific areas (the blue and red arches in the plot). This might 
lead to interesting correlations in rare decays of B s mesons. In fact the shape found in 
the right panel of figure 14 is familiar from the AM S constraint projected onto the plane 
(|(Mf 2 )Np|, arg(Mf 2 )Np), see e.g. figure 1 of [97]. This leads us to conclude that the NP 
amplitude in the present case is totally dominated by the tt contribution, and therefore 
governed by ^t(B s ). Indeed from the structure of in (9.5) it is easy to derive that 
the strong hierarchy in the matrix R(B q ), see appendix D, cannot be overcompensated by 
hierarchies in the right-handed quark mixing. 

10 A brief discussion of flavour observables 

This section is dedicated to the detailed discussion of flavour violating effects in the LRM. 
After analysing the possible size of NP effects in the various meson systems, we study the 
possible enhancements of the B s mixing phase and effects in Br (13 —> X q ^). Finally we 
turn our attention to the | V u b | problem, which was recently discussed in the context of RH 
currents [41-44]. 

10.1 A summary of possible size of NP effects in the different meson systems 

As discussed in section 9.3 the relative size of LR contributions in the various meson systems 
is given by the quantity ^t(q)- In figure 15 we show the possible sizes of |^i(g)| for the 
three meson systems as functions of each other for s = 0.1. We observe that in principle 
huge effects in the K system are possible, while the effects in the B s system are much 
more moderate and the effects in the Bj system are rather small. This indicates a rough 
hierarchy of NP effects Bj < B s <C K also expected for rare K and B^ s decays. Note that 
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Figure 15. The absolute values of ^t(g) for q — K,B^B S as functions of each other. The 
colour-code corresponds to the fine-tuning AgQ d . 
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Figure 16. The model independent correlation of A| L and (left panel) and the correlation of 
Ar s with 2(j) s (right panel) 

huge effects in the K system are possible only for large fine-tuning, since in that case the 
large tt contribution to K° — K° mixing has to be cancelled by other contributions. 

10.2 The phase of B s mixing and Br(S — > Xg'y) 

Now we want to briefly discuss the results of our analysis of observables related to the 
phase of B s mixing and Br(.B — > Xg'y). In the general case with and without fine-tuned 
points we observe no correlation between observables, except for model-independent ones 
such as A^-S^cf, and Ar s versus (fr s shown in figure 16. 

In figure 17 we show the correlation between Br(i? — > X s ^f) and Br(S — >• X^j), en- 
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Figure 17. The correlation of Br(B — > X s ^) and (Br(B — > Xdj)) showing the s dependence (left 
panel) and Re(V t ^) dependence (right panel) in colour. The red cross indicates the SM values, while 
the experimental central values are given by the dashed lines. 

coding the s dependence (left panel) and Re(V^) dependence (right panel, for s = 0.1) 
in colour. For s > 0.5 the values of Br(S — > X s ^f) cover a larger range that originates 
in the divergence of the function u(s) in the unphysical limit s — > l/\/2, see section 4.3. 
For small s the branching ratios depend only linearly on this parameter and the effects are 
much smaller. As can be seen from the explicit formulae given in section 4.6 Br(B — > X s j) 
depends linearly on the real part of V^. 

Combining the information of both plots in figure 17 we conclude that in order to 
enhance Br(i? — > X s j), we would need both s > 0.5 and a dominantly real and positive 
V^. As the SM prediction lies somewhat below the experimental value, albeit still in good 
agreement, a positive NP contribution to Br(B — > X s j) is welcome. 

Concerning the CP asymmetries in b — > we find moderate changes relative to the 
SM prediction 

A C p(b^djf u = -9.2%, (10.1) 
A C p(b^ s 7 ) SM = 0.4%. (10.2) 

Similar to the branching ratios also the CP asymmetries are s dependent. For s = 0.1 we 
find 

- 13.8% < A CP (b di) < -7.9% , (10.3) 

0.3% < A CP (b 87) < 0.6% , (10.4) 

while for s = 0.6 the possible range is slightly larger 

- 14.2% < A C p(b -> d-y) < -8.3% , (10.5) 
-0.18% < A CP (b 87) < 1.1% . (10.6) 

These effects will be difficult to disentangle from the SM contribution due to the large 
non-perturbative uncertainties present in these observables. 
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Figure 18. Maximal value for sce 2 |V^| as a function of s, for different allowed fine-tunings. A 
complete solution to the \V u b\ problem requires values for see 2 \V^\ above the dashed black line. 



Note that in figure 17 and in (10.3)-(10.6) we always used the maximal allowed Higgs 
mass (see section 11). For lower Higgs mass the constraints on the parameter space become 
much stronger, predicting a much more hierarchical structure for V^. These constraints 
however do not suppress Br(B — > X Sj( ri) since the latter decays depend only on the diag- 
onal element V^. Therefore decreasing the heavy Higgs mass M# enhances the effects in 
Bv(B^X s4l ). 

10.3 The | Kb | problem in the LRM 

As already discussed in section 5 in the LRM the true value of \Vh\ is given by the inclusive 
value of \V u b\ in (5.10), which remains unaffected by the presence of RH currents. On the 
other hand the corrections from RH currents to the vector and axial-vector couplings 
could in principle allow to explain the different values of \V u b\ found in exclusive decays 
and B + -4 t + v t as given in (5.11) and (5.12), respectively [41-44]. A solution to the 
\V U b\ problem can be provided if -see 2 1 l^^l is in the ballpark of 0.6 x 10~ 3 . Since we can 
only restrict see 2 < 10~ 3 , while \V]%\ is not constrained if arbitrary fine-tuning is allowed, 
studying the s dependence of the constraints is mandatory. 

Figure 18 illustrates the situation in the LRM. In general the black dashed line, nec- 
essary for a complete solution of the \V U \,\ problem, cannot be reached in the LRM. The 
tension increases significantly if we consider only points with small fine-tuning. While we 
cannot exclude that there are some models with RH currents that could solve this problem, 
from the point of view of the LRM considered here, this is not the case and any value of 
\V u b\ among (5.10)-(5.11) could be the true value. 
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Figure 19. The allowed range for the Higgs mass Mjj as a function of s for w 5.8 TeV and 
s R = 0.8. 

11 The role of the heavy Higgs mass 

Up to now in our analysis we have set Mh = 16/y/l — 2s 2 TeV, this ensures staying in the 
perturbative regime of the coupling 03 in the Higgs potential while suppressing the Higgs 
contributions in AF = 2. In this section we want to investigate the lower bound on the 
Higgs mass from our constraints. Therefore we now allow 03 to vary between 0.1 and 8, 
or equivalently the heavy Higgs mass between 2/Vl - 2s 2 TeV and 16/Vl - 2s 2 TeV. The 
result can be seen in figure 19, in particular we find the 'soft' lower limit on Mjj to be 

2.4TeV^Mtf. (11.1) 

This limit is 'soft' in the sense that lower values are not rigorously forbidden, but our 
parameter scan did not reveal any points which allow for lower masses. One might suspect 
a very large fine-tuning for this low Higgs mass, but while the fine-tuning does increase 
with decreasing Higgs mass there are still points with fine-tuning AgQ d < 10 even for the 
lowest possible mass. This is of course also a test of the idea of fine-tuning. Finally for low 
Higgs masses the matrix V R shows a very hierarchical structure as can be expected. The 
mixing pattern is roughly given by 

s 12 ~ ©(lO" 2 ) , si 3 ~0(HT 4 ), s 2 3~0(l(T 3 ). (11.2) 
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12 Comparison to other models 



12.1 Preliminaries 

Now we compare the results found in the LRM to those obtained in other NP models. 
The first part of this section is dedicated to the comparison with two specific versions of 
LR models discussed in the literature. Subsequently we compare the pattern of flavour 
violation in the general LRM with patterns found by us in a number of extensions of the 
SM. A review of those analyses can be found in [6]. Our discussion below deals only with 
AF = 2 transitions and B — > decays. 

In order to increase the transparency of the comparison in question let us summarize 
the most characteristic features of NP contributions in the LRM: 

• The dominance of the operator in all AF = 2 transitions: K° — K°, B® — B® 
and B® — B®, with the biggest impact on the K° — K° system. 

• While both box diagrams with Wl — Wr exchanges and tree level neutral Higgs 
contributions are the sources of the effects in AF = 2 observables, the latter contri- 
butions are most important. Governed by the same operators and having the same 
quark mixing factors the trees always dominate over the boxes. 

• The contributions of induced RH couplings of Wl and of to the B — > X s ^ 
decays are enhanced by the ratio mtjm^ due to the chirality flip on the top quark 
propagator. Moreover contributions receive an additional enhancement for large 
s. 

• All effects listed above are beyond the usual MFV framework. 

With this picture in mind we are in the position to make a transparent comparison 
with other versions of LR models and with a number of extensions of the SM. 

12.2 LR models with additional flavour symmetries 

Recently two different approaches have been presented to reconcile the LR models with 
flavour data even for relatively small symmetry breaking scales. In both cases this is 
achieved by introducing an additional symmetry that removes all flavour changing Higgs 
effects from the down quark sector. 

In a first paper Guadagnoli and Mohapatra [158] suggested to extend the Higgs sector 
by a discrete Z4 symmetry. Eventually this leads to removing all flavour violating Higgs 
couplings from the down quark sector and the strongest constraints on the heavy Higgs 
mass then comes from D° — D° mixing. AS = 2 and AB = 2 processes are then governed 
by Wl — Wr box diagrams. While in [158] the manifest version of the LRSM has been 
considered, it is straightforward to extend this model to the general LRM leading to a 
potentially interesting phenomenology. 

A different approach to avoid the stringent flavour constraints has been proposed in 
[62]. Here a gauged LR symmetric flavour group SU(3)q l x SU(3)q r has been introduced 
following the idea of Grinstein et al. [82] . In this scenario due to a different structure of the 
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Higgs sector only the SM Higgs boson is present and no dangerous tree level FCNH effects 
arise. The potentially dangerous tree level contributions from flavour gauge bosons are nat- 
urally suppressed by a see-saw-like mechanism, so that the most relevant new contributions 
to AF = 2 processes arise again from Wl — Wr box diagrams. 

The effects in the B — > X S( fY decays turn out to be small in this scenario. This is 
related to both the absence of a heavy charged Higgs and the fact that in this model Wl 
and Wr do not mix with each other. Consequently no chiral enhancement cx mt/rrib is 
possible in this case as in the scenario considered by us. On the other hand as pointed 
out in [102] in such models new flavour violating neutral gauge bosons could provide an 
enhancement of B — > X s dJ through a chiral enhancement cx Trip/mi,, where F denotes a 
new vectorial neutral heavy fermion. Yet, as shown in [102] in models in which fermion 
masses are generated by a see-saw mechanism such contributions are strongly suppressed 
by the mixing between SM quarks and heavy fermions so that finally the NP contributions 
remain small. 

12.3 2HDM with flavour blind phases: 2HDM MfV 

Concerning the AF = 2 transitions the most spectacular difference in 2HDM MFV [75, 76] 
from the LRM considered here is the absence of relevant direct contributions to ek- The 
reason is that in 2HDM MFV the flavour violating neutral Higgs couplings are effectively 
generated at one-loop level and are necessarily proportional to the masses of external quarks 
(md, s in this case) and not to up-quark masses as in the case of the LRM, where flavour 
violating Higgs couplings are present in the fundamental Lagrangian. 

Therefore in 2HDM MFV can only be found close to the data if sin 2/3 ~ 0.80 which 
corresponds to scenario 2 in section 7. The interplay of flavour blind phases with the CKM 
phase allows then to bring S^Ks t° As experimental value, while enhancing automatically 

to values close to 0.3. Having less parameters than the LRM, this prediction is rather 
unique. Finding experimentally to be negative would rule out 2HDM MFV . For LRM 
this would still not be a problem as one would move then to a different oasis in the parameter 
space of this model than considered by us in 7.3.2. 

It should also be emphasized that the neutral Higgs masses in the 2HDM MFV are of 
the order of a few hundred GeV, whereas in the LRM they are by at least one order of 
magnitude larger. Therefore LHC should easily distinguish between these two scenarios. 

12.4 Right-handed MFV 

Recently, in [44] , we have studied the effects of RH currents by means of an effective theory 
approach under incorporation of an extended MFV principle. 

The similar symmetry pattern with respect to the LRM could lead to the premature 
assumption that both models can be matched by integrating out the heavy particles in 
the LRM. However as we have already mentioned in section 3.8, the operator structure 
of both models turns out to be complementary. Furthermore, due to the extended MFV 
mechanism the specific flavour structure precludes this possibility. 

The differences display also in the phenomenology of both models. Even though the 
general parametrization for the RH mixing matrix is given by the same parameters, after 
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imposing all constraints the matrices display a different pattern. This can be seen explicitly 
when comparing figure 12 to the structure of the matrix given in equation (105) of [44]. 
The fact that the \ V u b\ problem cannot be solved in the LRM whereas this was the case in 
RHMFV demonstrates once more the difference in the RH mixing matrices. 

Let us have a closer look at further flavour observables. While the RHMFV was 
originally designed to explain a large S^Si a full theory like the LRM allows a more flexible 
analysis. Here no particular assumptions about S^^ have been made and values within 
the full range are allowed. Since in the LRM we take into account the variation of all 
phases, the high number of parameters allows to find regions in parameters space where all 
SM tensions can be resolved. This was not the case in RHMFV where the S^k s -sin(2/3) 
tension persists under the assumption that S^hcb is significantly enhanced over its SM value. 

On the other hand the general pattern of NP effects within the different meson sectors 
show similarities for both models. For example in both cases we find the rough hierarchy 
B<i < B s <C K for effects in the corresponding observables. The observable ek acts in both 
cases as strongest constraint. Finally both frameworks cannot explain the anomalous Zbb 
coupling. 

12.5 Randall-Sundrum with custodial protection 

In the RS model with custodial protection (RSc), AF = 2 operators are generated already 
at the tree level by the exchange of heavy Kaluza-Klein gluon and electroweak gauge boson 
modes (see [68] for an extensive phenomenological analysis). Tree level flavour changing 
couplings of the SM Higgs generate AF = 2 transitions as well, but they turn out to be 
subleading with respect to the aforementioned KK gauge boson contributions [159, 160]. 
In contrast to the LRM studied here, only K — K mixing is dominated by the chirally 
enhanced operator generated by KK gluon exchange. Since the chiral enhancement is 
absent in the case of Bd, s — Bd,s mixing the latter processes are dominated by the operator 
QY LL and electroweak KK modes are relevant. Consequently while the constraint from 
ek generically puts strong bounds on the model in question [155], the contributions to 
M 1 2* are suppressed by the RS-GIM mechanism [161]. Still non-zero effects are generally 
expected, so that a solution of the ek - S^Ks tension is possible within this framework. 
Also can receive large enhancements. 

In contrast to AF = 2 transitions, the dipole operators governing the B — > X s dj 
decays are generated first at the one loop level. As in the LRM the mj, suppression can 
be overcome by a chirality flip inside the loop. However in contrast to that model, in the 
RSc the primed operators C 7 SG are generally dominant [161-163], so that an enhanced 
branching ratio is generally expected. 

12.6 Four generations: SM4 

As opposed to the scenarios discussed so far, there are only SM operators contributing to 
AF = 2 processes and the NP effects are primarily enhanced through the non-decoupling 
effects of the i'-quark. The presence of new mixing angles and of new phases allows then 
to solve the S^k s ~ £ K anomaly in a straightforward manner and simultaneously enhance 
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Sfc/, if necessary. Moreover the agreement of theory and data can be improved in various 
observables. 

Concerning the B —> X s j decay, as seen in figure 17 of [69], the branching ratio can be 
enhanced up to 4- 10 -4 but also significantly suppressed below 3- 10 -4 . Large enhancements 
like the ones from in the LRM are not possible in the SM4, in particular if one wants 
to get interesting deviations from the SM3 predictions for and B s — > ^ ~fj,~ . 

12.7 Littlest Higgs with T-parity 

As in the case of the SM4 no new operators enter in the LHT model and FCNC processes are 
generated first at the one- loop level [67, 164]. Still due to the absence of an intrinsic flavour 
structure in the mirror fermion sector the effects on AF = 2 transitions are generally large, 
and in particular ek puts strong constraints on the new mixing parameters. At the same 
time the presence of new mixings and phases in the mirror sector allows to solve anomalies 
including a possible enhancement of S^. However the effects are not as strong as in the 
case of SM4 [165, 166]. Due to the absence of right-handed flavour violating currents, no 
chirality enhanced contribution to B — > X s ^f arises and the effects remain very small. 

13 Summary 

In this paper we have presented a complete study of AS = 2 and AB = 2 processes in 
the LRM. This includes e K , AM K , AM S , AM d , A q SL , AF q , and the mixing induced CP 
asymmetries S^k s and S^. We have included the new contributions from box diagrams 
with Wr gauge boson and charged Higgs (H^) exchanges and tree level contributions from 
heavy neutral Higgs exchanges. We have also analysed the B — > Xg^ decays that receive 
important new contributions from the Wl — Wr mixing and exchanges. Compared to 
the existing literature the novel feature of our analysis is the search for correlations between 
various observables, simultaneous inclusion of all relevant contributions, in particular Higgs 
contributions and an improved treatment of QCD corrections. Our main findings are as 
follows: 

• We find that the LRM is put under pressure by the ek constraint. This is due to 
the tree-level contributions of the neutral Higgs scalars and the related LR oper- 
ators whose contributions are enhanced through renormalisation group effects and 
their chirally enhanced hadronic matrix elements. While this problem has been 
known for many years, we stress that even if the heavy Higgs masses are of order 
15 — 20 TeV, these contributions are not only important but even dominant. In- 
creasing these masses much more would make the Higgs system non-perturbative. 
Leaving them out from phenomenological analysis, as done in some papers, is simply 
wrong. Figure 3 makes this point explicitly. On the other hand we show that there 
are structures (hierarchical or inverted) of the matrix where the ex-constraint 
can be satisfied without large fine-tuning of parameters. 

• The contributions of Wr to AF = 2 observables, in particular ek, is important but 
in view of the new improved lower experimental bound on its mass much less prob- 
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lematic than the tree level Higgs contributions, although also here the LR operators 
mentioned above enter. Consequently, if only these contributions were present, the 
matrix V R would be somewhat less hierarchical than the CKM matrix. 

• The charged Higgs being as heavy as the neutral ones plays a subdominant role 
in AF = 2 processes but can give significant contributions to B — > X Si dl decays. 
These contributions are often neglected. We stress that this is not justified. Figure 4 
makes this point explicitly. We find that the branching ratios for these decays can be 
significantly affected by new contributions but the CP-violating effects are too small 
to be distinguished from the SM results in view of hadronic uncertainties. 

• In our analysis of FCNC processes we have taken into account all existing constraints 
from electroweak precision observables and tree level decays. In this context we have 
found that in these models the SM problem with RH Z — )• bb couplings and the 
so-called \V u b\ -problem cannot be solved. The increased value of M\y R due to LHC 
combined with FCNC constraints does not allow to explain the difference in exclusive 
and inclusive determinations of \V u b\ with the help of RH currents within the context 
of the LRM considered here. 

• We have found that the NP effects are largest in the K system but large effects can 
also be found in the B s system. Significant but smaller effects are found in the Bj 
system. 

• Guided by the persisting S^k s ~ £ K tension present in the SM we have constructed a 
simple analytic expression for V R , which is given in (7.5). It is given in terms of two 
mixing angles S13 and S23 and two corresponding phases <f>\ = <fii3 and 02 = 023 arid 
allows to remove or soften the tension in question for scales M\y R ~ 2 — 3 TeV in the 
reach of the LHC. In this scenario for V R interesting correlations with CP violation 
in the B s system are present that depend on the value of \ V u b\ chosen as well as the 
numerical values of the elements of V R . In two \V u b\ scenarios considered by us 
can be significantly enhanced over the SM value. When the data on will improve 
one will be able to learn more about the structure of V R . 

• A more involved numerical analysis allows to find other structures of V R with all its 
nine parameters entering the game that without large fine-tuning of these parameters 
allows to obtain interesting results. In particular as shown in figure 12 three oases are 
found in the space (|V^|, |V^|, |V C R |), of which two have low fine-tuning. The pattern 
of flavour violation in the scenario with "normal" hierarchy in has been presented 
in several plots in section 9 and the corresponding Wolfenstein-like parametrization 
has been derived. 

• We find a soft lower limit for the mass of the heavy Higgses of Mh ^ 2.4 TeV, see 
figure 19. 

Our analysis of the LRM presented in this paper was dominated by the observables 
on which already good data exist. We have seen that already these data put stringent 
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constraints on this NP scenario. However there are other observables not discussed by 
us, that have not been measured yet but the experimental upper bounds on them put 
significant constraints on the LRM. Among them a prominent role is played by the neutron 
electric dipole moment (EDM). The most recent analysis in [28] demonstrates that the 
neutron EDM provides lower bounds on W and heavy neutral Higgs masses that are 
competitive with those coming from Ek- Determining what impact the present bounds on 
the neutron EDM have on the LRM considered by us would require a detailed analysis of 
all contributions and taking into account the significant hadronic uncertainties. We leave 
such an analysis for future work. Similarly, in this decade other observables, like rare B 
and K decay branching ratios and related CP violating observables, will be measured with 
sufficient precision so that we will be able to find out whether this NP scenario is viable. 
Of course the discovery of new gauge bosons W and Z' at the LHC would be the most 
spectacular manifestation of the kind of new physics analysed in our paper. 
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A Higgs potential 

The most general renormalisable Higgs potential invariant under parity is given by [28, 29, 
167, 168] 
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where there are a total of 18 parameters, Hi 2 31 -^1,2,3,4, Pi, 2,3,4, 0^1,2,3, and /?i,2,3- Due to 
the left-right symmetry, only one of them, the coupling involving 02 can become complex 
and all other couplings are real. 

We have investigated the possibility of a more general Higgs potential, where no left- 
right symmetry is imposed. However as the parameters of the Higgs potential do not 
affect our flavour analysis, taking into account that leading order couplings of Higgs and 
quark fields are independent of these parameters, we find it sufficient to give here the left- 
right symmetric form. In fact among the parameters of the Higgs sector only 03 in (A.l), 
determining the heavy Higgs masses, is relevant for our analysis. 

B Gauge boson masses and mixings 

The gauge boson mass matrices can straightforwardly be obtained from the relevant terms 
of equations (2.16) and (2.17). Diagonalising the resulting mass matrix for the charged 
gauge bosons yields the mass eigenstates 

W ± = Wf + sce^ aS ^e 2 W± , (B.l) 
s w 

W* =W£- sce ±iaS ^e 2 W± , (B.2) 
sw 

where we defined v = y k 2 + k' 2 , c = k/v, s = k'/v and we neglected higher order terms 
in e = v/kr. We also introduced the mixing angles 



9 12 S R -1 2 

sr = j , c R = \l 1 - s R , s w = , c w = W 1 - sfc . 



9' 2 + g 2 R " J{g L /g R ) 2 + s 2 R 



(B.3) 



The corresponding masses read 



e 2 v 2 



(Mw? = (1 - 2s 2 c 2 e 2 ) , (B.4) 



^s w 

„2,J2 



C W S R 



1 



{M Wr Y = ^{l + -e 2 \ . (B.5) 



For the neutral gauge bosons we can write 

A = s w Wl + s r c w Wr + c R c w B , (B.6) 

Z = c w W[ - srs w (l - £§- e 2 ) W R - crsw (l + J^e 2 ) B (B.7) 

Z> = - S -^K e 2 Wi + c R (l + ^ej W R -s R (l- ^e 2 ) B . (B.8) 



The masses are given by 



(M A ) 2 = 0, (B.9) 

(Mz)2 = ^{ i - r ¥) ■ < rio » 

^ b W c W V ^ / 



s R c r c w V * 
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C Goldstone boson and Higgs mass eigenstates 



The doubly charged components <5^jj lead to two physical doubly charged Higgses. 

The singly charged fields cj)^ 2 and 5^ R form the Goldstone bosons of and W 1 ^ 1 
and two charged Higgses. The Goldstone bosons are given by 



G 
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G 



±i 



c 1 



A 2 ) 



- se Tia (1 


-cV) 


se ±ia , 
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V2 ^ 


;= e 



V2cse^ Q e<5 



The physical singly charged Higgs states are linear combinations of 5^ and 



se 



„2\2 



t>f+C 



(c 2 



„2\2 



V2 



-eS 



Ft 



Finally the neutral Goldstone boson and Higgs fields are built out of /A ° 



1,2 



Defining tt° = clm<^ — slm(e 



~ %a 4>2), the Goldstone bosons of Z and Z' read 



(CI) 
(C.2) 
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and 8 L R . 
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C -R c 2 
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lmS° R 
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(C.4) 
(C.5) 



In addition there are six neutral Higgs fields in the spectrum, that are linear combinations 
of slm^ + clm(e- ia (/) 2 ), Re0?, Re(e- ia <$), Re<5^, ReS° L and Im<5£. 

In order to make a more detailed statement about the Higgs mass eigenstates one has 
to diagonalize the Higgs potential (A.l). Considering the neutral 8 by 8 mass matrix, its 
diagonalisation can be best done by using perturbation theory. While in [28] a further 
hierarchy, namely k' <C k, has been assumed, we do not restrict ourselves to this case but 
determine the Higgs mass eigenstates to leading order for arbitrary s = k' /v. It turns 
out that the leading order Higgs couplings are not sensitive to the detailed structure of 
the potential and in particular the parity invariant potential yields the same result. The 
leading order mass eigenstates of the Higgs fields of our interest are then given by 



h° = V2 (cRe0? + sRe(e- ia ^)) 
Hi = V2 (sRe^ + cRe(e" ia ^)) 
ff° = ^2 (slm^l + clm{e- ia ctP 2 )) 



(C.6) 
(C.7) 
(C.8) 



where ho can be identified as the light Higgs with a mass of 0(v), while Hf and H® are two 
new flavour-violating neutral Higgses with masses O(kr). To leading order their masses 
are equal to each other 



2 

M 2 _ M 2 _ M 2 _ a 3 K R 



a 3 



(C.9) 



In the course of the numerical analysis we regard Mh as a free parameter. These results 
agree with the ones given by the authors of [28] in the limit of small s<l. A more explicit 
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analysis of the Higgs sector of this class of models can be also found in [28, 167, 168]. 
Charged Higgs effects in LR models are often neglected in the literature with the argument 
that they have to be small in FCNC processes as they have to take place at one loop order 
[27]. This is in fact true for the AF = 2 processes but not for B — > X s ^- This has already 
been pointed out in [54, 55, 58-61] in LR models and confirmed here by us. As in the case 
of the neutral Higgs sector we are only interested in leading order couplings and masses. 
We find that the mass of the lightest charged Higgs, 

H± = se ±ia (j)f + c^f , (CIO) 

is given to a very good approximation by 

M H+ = M H (C.ll) 

and hence cannot be chosen to be as light as 1 TeV as done sometimes in the literature. 



D Numerical details for AF = 2 



Here we give the values of Rij as defined in (3.40) and used in the main part of our numerical 
analysis. Including all contributions we find for e = 0.03, sr = 0.8 and s = 0.1 
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For s = 0.5 the Higgs contributions get enhanced and the calculation of the Rij yields 
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Keeping only gauge boson contributions we find (no s dependence) 
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The contributions are significantly smaller than gauge contributions. 
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E Feynman rules 



In this part of the appendix we present the Feynman rules up to and including 0(e 2 ) 
corrections, except for the couplings to Higgs bosons where we restrict ourselves to the 
leading order couplings 10 . The Feynman rules are given in mass eigenstates of the particular 
fields. We want to stress that we keep here the correct notation W and W', as introduced 
in appendix B, in spite of using throughout the phenomenological analysis the simplified 
notation Wl and Wr for the mass eigenstates. 

Concerning triple gauge couplings a comment is in order. The Dirac structure of all 
vertices is the same, 




C [v»v( k - V)p + Vupip - q)n + VpniQ - k )u] - 



where Vl~ = W+, W'+, V u = W u , W' v , V® = A p , Z p , , Z' and k,p, q are their incoming 
momenta. Therefore in table 12 we collect only the coefficients C of the respective couplings. 
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Table 7. Fermion couplings to the massless gauge bosons X^: photon and gluons G" 
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Table 8. Fermion couplings to W + and W' + 



3 We would like to thank Jennifer Girrbach for checking the Feynman rules presented here. 
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Table 9. Fermion couplings to neutral gauge bosons Z and Z 
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Table 10. Fermion couplings to charged Goldstone bosons 
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Table 11. Fermion couplings to neutral Goldstone bosons G° and G'° 
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Table 12. Triple gauge couplings, involving a Z and Z' boson and the photon, respectively 
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Table 13. Charged Goldstone couplings to the photon and the Z boson 
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Table 14. Charged Goldstone couplings to the Z' boson 
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Table 15. Couplings of charged and neutral Goldstone Bosons to the W and W boson 
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Table 16. Couplings of charged gauge bosons to the G° and G'° boson 
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Table 17. Couplings of charged and neutral gauge bosons to the G + and G' + boson 
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Table 18. Fermion couplings to the flavour violating neutral Higgses. Here mj and m d denote the 
ath up and down quark mass, respectively. Summation over a is understood. 



v? L d R H+ 




i 


( m l V- - 


- 2cse~ ia V}-m>^\ 


(i 


-2s 2 )v 


u R d j L H + 






V£m J d - 


2cse- ia m i u V^j 


(i- 


2s 2 )v ( 



Table 19. Charged Higgs couplings to fermions. Here m£ and m d denote the ath up and down 
quark mass, respectively. 
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